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PREFACIO A LA EDICIÓN EN LENGUA ESPAÑOLA 


A 


Con motivo de la decisión de publicar esta obra en lengua 
española desearía dirigirme al lector con algunas palabras. 

Hace ya cerca de veinte años que apareció este libro. En aquel 
entonces la programación lineal se encontraba en el proceso 
de su desarrollo, estábamos apasionados con las cuestiones de su 
teoría, aplicación y divulgación. En este pequeño libro traté 
de explicar los conceptos básicos de la asignatura y la esencia 
de algunos de sus métodos de cómputo. 

El libro se agotó rápidamente, continuaron las investigaciones 
sobre el tema; han aparecido muchas obras en las cuales, bajo 
mi opinión, se explica mejor parte del material. A pesar de eso, 
hace relativamente poco tiempo, este libro fue reeditado en Japón, 
en la RDA, ha sido traducido a las lenguas de otros paises. 
Ahora la editorial “MIR” está preparando en lengua española 
la edición de la serie completa de libros “Lecciones de divulgación 
de matemáticas” en la que éste salió. La traducción de este libro 
se ha hecho integramente por dicha edición. Observaremos que 
desde aquel tiempo las computadoras se han perfeccionado 
muchísimo y por eso las referencias a los ejemplos de rapidez 
de resolución de algunos problemas, así como las dimensiones 
de éstos, deben ser interpretados por el lector como factores que 
caracterizaron el período inicial de la aplicación de la programación 
lineal. Actualmente las dimensiones de los problemas que se tratan 
son mucho mayores y las velocidades de su resolución mucho más 
altas. Las explicaciones que se dan sobre los conceptos y definiciones 
fundamentales de la programación lineal y las bases de los métodos 
de cálculo por el criterio de costo y de tiempo, pueden hoy 
en día abrir al lector el camino a esta interesante asignatura 
de la matemática moderna, a la programación lineal. 

Estaré satisfecho de que el lector encuentre útil para sí este 
pequeño libro y al mismo tiempo le expreso mi simpatía. 


5 de mayo de 1976 A. S. Bársov 


PREFACIO 


En este libro se examinan cuestiones de la teoría y de los 
métodos de resolución de algunos problemas de programación 
lineal. El está destinado a un amplio círculo de personas ocupadas 
en el empleo de métodos matemáticos en la organización 
y la planificación de la industria. 

Se estudian los fundamentos de la programación lineal. 
Al hacerlo se presentan tan sólo los datos y demostraciones que 
son necesarios para una exposición elemental de los métodos 
de programación lineal. 

El trabajo se realizó a base de las conferencias dadas por 
el autor en el año 1957 para personas que se dedicaban 
a la resolución de problemas de programación lineal en las 
máquinas computadoras electrónicas. 

L. A. Lyusternik, miembro correspondiente de la A.C. de la 
U.R.S.S., en el año 1959 revisó con atención el material de las 
conferencias, dio una serie de valiosos consejos y contribuyó 
a la edición del presente trabajo. 

El autor agradece a los profesores A. A. Liapunov y N.S. Kra- 
silnikov por su ayuda en la solución de las dificultades que 
surgieron en el proceso de la preparación de este libro. 

El autor está particularmente agradecido al redactor V. D. Ro- 
senko por su minuciosa labor que sirvió considerablemente para 
mejorar el libro. 


A. S. Bársov 
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INTRODUCCIÓN 


La tarea del desarrollo ulterior de las fuerzas productivas, 
de la mejora de la planificación de la industria socialista y el 
aumento de la efectividad económica de las inversiones básicas 
en nuestro país adquiere de año en año cada vez mayor 
importancia. 

La diversidad de posibles soluciones técnicas y caminos 
de desarrollo en la industria actual, las interrelaciones entre las 
diferentes ramas de la economía nacional y otros problemas 
económicos hacen que las tareas planteadas anteriormente sean 
excepcionalmente difíciles. 

Para solucionar estos problemas, los métodos matemáticos 
y en particular el método de programación lineal, así como los 
medios técnicos modernos, las máquinas computadoras electrónicas 
pueden prestar una ayuda sustancial. 

La teoría de la programación lineal, que surgió hace dos 
décadas, actualmente ha obtenido una amplia utilización práctica 
particularmente en el terreno de la organización y planificación de 
la industria. 

Los primeros trabajos en este sentido fueron los de L. V. Kan- 
torovich, miembro correspondiente de la A.C. de la U.R.S.S.*. 
En esos trabajos se expusieron métodos matemáticos para resolver 
problemas tales como el del aumento de la efectividad del 
transporte, el cálculo de los regímenes óptimos de producción, 
la distribución racional de los materiales industriales, etc. 

Posteriormente fueron creados métodos generales de pro- 
gramación lineal tales como, por ejemplo, el simplex, el combina- 
torio y otros métodos que se emplean eficazmente en la solución 
de diversos problemas para determinar el óptimo en la planifi- 
cación. Dantzing, Charnes y una serie de cientificos soviéticos 
y extranjeros se ocuparon de la elaboración de estos métodos. 

La programación lineal abarca métodos de solución de problemas 
de óptimo en los que hay muchas variables relacionadas entre 
si y subordinadas a unas determinadas condiciones de restricción. 
El planteamiento de los problemas de programación lineal se puede 
formular de la siguiente manera: 

Tenemos cierta magnitud (por ejemplo el costo, el tiempo) 


*) L. V. Kantorovich hoy es miembro efectivo de la 
A. C. de la U.R.S.S. y laureado con el Premio Nobel. 
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que es función lineal de una serie de variables. A su vez, las 
variables tienen que satisfacer a las restricciones expresadas en 
forma de un sistema de desigualdades o igualdades lineales. 

Hay que buscar aquellos valores no negativos de las variables 
con los que la magnitud que sea su función adquiera el valor 
menor (el mayor). 

En calidad de ejemplo veamos el problema de transporte. 
Este problema se puede formular de la forma siguiente: 

De m puntos de partida dados, en cada uno de los cuales 
hay a; unidades de cargamento, se deben llevar b; unidades 
de la carga indicada a cada uno de los n puntos de llegada 


A A EA 


Se exige planificar el transporte de tal modo que los gastos 
sean mínimos. Designaremos con x;; la cantidad de carga trans- 
portada del punto de partida ¡ al punto de llegada j. En tal 
caso el contenido matemático del problema se reduce a encontrar 
los valores no negativos de x,, que satisfagan las ecuaciones 


y con los cuales el costo general del transporte 


sea el menor. Aquí c;,; es el costo de transporte de una unidad 
de carga desde el punto de partida ¡ hasta el punto de llegada j. 

En un serie de casos prácticos importantes, el problema se 
plantea así: se necesita planificar el transporte de cargamentos 
desde m puntos de partida dados hasta n puntos de llegada 
de tal forma que todas las operaciones de transporte se realicen 
en un plazo mínimo. 

- Otro ejemplo del empleo de la programación lineal puede ser 
el problema siguiente: 

En muchas fábricas, la producción de diferentes artículos se 
realiza en líneas automáticas. En estos casos pueden surgir diferentes 
problemas referentes a la organización más racional de la producción. 

Supongamos, por ejemplo, que un taller dispone de m máquinas 
para fabricar n artículos distintos. Cada máquina i (i = 1, 2, ..., m) 
se caracteriza por el tiempo posible de trabajo mensual b,, 
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la norma de tiempo f;; para elaborar una unidad del artículo j 
(¡=1, 2, .., n) y por el costo c;; de fabricación de una unidad 
del artículo ¡ en la máquina i. Si al taller se le ha dado la tarea 
de sacar en el mes próximo una determinada cantidad a, de cada 
uno de los diferentes artículos, surge el problema de organizar 
el trabajo en tal forma que esta tarea se cumpla con los gastos 
mínimos. Si designamos con x;; la cantidad de artículos j que 
se fabrican en la máquina i, el problema planteado se reducirá 
a una distribución del trabajo entre las máquinas que satisfaga 
las condiciones 


n 
yA Xij 2 M5 > ¡ij < by 
¡=1 j=1 


y que lleve el valor total de los costos 


al mínimo posible. 

En los problemas de programación lineal las condiciones que 
se aplican a la zona de valores admisibles de las variables se 
determinan por un sistema de desigualdades o igualdades lineales. 
Al mismo tiempo, la función cuyo valor mínimo (máximo) se busca, 
es también una función lineal de las mismas variables. Esto lo 
subraya la propia denominación de programación lineal. 

Para la determinación de la solución óptima, los métodos 
de programación lineal exigen el examen de varias soluciones. 
Al analizar problemas prácticos, por ejemplo, el problema del 
empleo racional de la maquinaria o de las empresas con determinadas 
condiciones restringentes, el paso de una solución a otra corresponde 
a un estudio sucesivo de diferentes programas de producción. 
De aquí viene el nombre de programación lineal. 

El problema de la programación lineal es un problema 
de búsqueda de un punto de cierta región en el que la función 
adquiere el valor máximo (mínimo). Por eso surge la pregunta 
natural: ¿Por qué no es suficiente en este caso emplear los 
conocidos métodos clásicos de resolución de los problemas de 
cálculo de valores extremos de las funciones, por ejemplo, el 
método de Lagrange? 

La causa reside en que los métodos clásicos exigen la existencia 
de derivadas parciales de la función en el punto en el que se 
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alcanza el extremo mientras que la función lineal llega a su valor 
extremo en el límite de la región, donde las derivadas parciales 
no existen. 

Eso fue lo que sirvió de motivo para la creación de nuevos 
métodos de resolución de problemas de extremo entre los que se 
encuentra la programación lineal. 

La práctica de resolución de problemas de programación lineal 
muestra que cuando es grande el número de variables, para 
resolver tales problemas es necesario emplear computadoras 
electrónicas. La máquina resuelve en unos dos o tres minutos 
problemas en los que el hombre gastaría hasta una semana. 
Al ser muy grande el número de las variables, estos problemas 
pueden ser resueltos solamente por medio de computadoras 
electrónicas. 

Puede servir de ejemplo el cálculo del plan óptimo de transporte 
de arena de construcción a las obras de edificación de la ciudad 
de Moscú. En este problema había 10 puntos de partida y 230 puntos 
de llegada. El plan óptimo de transporte calculado en la computadora 
“Strelá”, en diez días del mes de junio de 1958 dio un 11% 
de economía. 

A continuación se examinan los fundamentos matemáticos y los 
métodos de resolución de algunos problemas de programación 
lineal y, en particular, de los problemas de transporte. 
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CAPÍTULO 1 


ALGUNOS CONCEPTOS Y DEFINICIONES 
DEL ALGEBRA LINEAL 


En este capítulo se exponen los conceptos básicos y las 
definiciones del álgebra lineal de espacio m-dimensional que son 
necesarios para la resolución de problemas de programación lineal. 


$ 1. CONCEPTO DE ESPACIO M-DIMENSIONAL 


Cualquier trío ordenado (a,, a,, az) de números reales puede 
ser interpretado geométricamente como un punto en el espacio. 
En concordancia con esta representación geométrica se ha tomado 
en la matemática la definición siguiente: el espacio tridimensional 
es el conjunto de todos los posibles trios ordenados (a,, a,, ax) 
de números reales *). También se dice que el sistema de números 
(a,, a,, az) determina el punto M en un espacio tridimensional 
con las coordenadas a,, a,, az o el vector P con los componentes 
4;, Az, Az en ese mismo espacio. Estos tres números reales son 
insuficientes para representar ciertos objetos, procesos o estados. 
Por ejemplo, para determinar la posición de un cuerpo sólido 
en el espacio son necesarias seis coordenadas. En el caso de que 
en una zona económica se produzcan determinados artículos 
industriales y agrícolas, como vagones de ferrocarril, automóviles, 
trigo, leche, cerillas y otros, entonces para la característica de esa 
producción industrial y agrícola de la zona se necesitará una 
sucesión ordenada de números reales. Por ejemplo, la tabla 1 
indica que la zona Ne 2 produce anualmente a,, toneladas de carbón, 
a,, toneladas de mineral de hierro, a,3 toneladas de acero, ..., 4», 
toneladas de trigo. 

Similarmente se definirán, también por medio de sucesiones 
ordenadas de números, por ejemplo: la cantidad de combustible 
de aviación de diferentes clases que se emplea en tal pais, 


*) Por analogía el conjunto de todos los posibles 
números reales (a,) es un espacio unidimensional cuyo modelo geométrico 
puede ser la recta; el conjunto de todos los posibles pares de números 
reales (a,, a,) es un espacio bidimensional cuyo modelo geométrico puede 
ser la superficie plana. 
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la cantidad de artículos de diferentes clases que se encuentran 
en tal almacén, etc. 
Tabla 1 


Carbón | Miner. 
d 
(t) 


Zona N2 1 


Estos ejemplos muestran la utilidad del estudio de todos los 
posibles conjuntos de sucesiones ordenadas de m números reales 
para el caso de que m sea cualquier número entero positivo. 
Como se sabe, un sistema ordenado de m números reales 


(41, Az, ..., Aj, «.., Am) se lama vector m-dimensional. Los números a,, 
¡ = 1, 2,..., m, se denominan componentes del vector P (a,, a,, ..., 4). 
Los vectores P (a,, az, ..., Am) y Q(b,, b», ..., b,y) se consideran 


iguales si y sólo si coinciden sus componentes que se encuentran 
en los mismos lugares, o sea si a; = b; para todos los ¡ = 1, 2, ...,.m. 

En el caso de que nos interesara conocer la productividad 
total de las diversas clases de producción de dos zonas económicas, 
evidentemente ésta podría calcularse sumando las productividades 
respectivas de estas zonas.*? Si la productividad de la zona Ne 1 
en todas las clases de productos se determina por medio del vector 
P, (a, ;, 412, ..., 41m) y la de la zona Ne 2, por el vector P, (a,;, a», 
..., 7), entonces la productividad total de estas dos zonas estará 
caracterizada por el vector 


OQ (411 + 421, 412 + 022, «> Aim + 02m): 


S1 la productividad de una zona se determina por el vector 
P =P(a,, a, ..., Am), al incrementarse la productividad de cada 
uno de los productos en k veces, el aumento de la productividad 
de la zona se podrá expresar con el vector Q, el que será 


Q =Q (ka,, ka», ..., KA). 


*) Comprenderemos en este caso por productividad 
la cantidad de producto obtenido en un determinado lapso. 
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Las definiciones introducidas son sencillamente una generalización 
de las operaciones ya conocidas con los vectores en un espacio 
tridimensional. Al aplicar el concepto del espacio tridimensional 
a la sucesión de números reales m, obtendremos la importante 
definición siguiente: el conjunto de todos los posibles m-dimensionales 
vectores P (a,, a», ..., 4, ) con componentes reales se llama espacio 
m-dimensional y se designa P”. Según la definición, sumar dos 
vectores m-dimensionales consiste en obtener un tercer vector 
R=P+0Q con sus componentes iguales a las sumas de los 
respectivos componentes de los vectores sumados. Multiplicar el 
vector P por el número k quiere decir multiplicar cada componente 
por ese número. 

Se dice que el vector P (a,, a,, ..., a,) es proporcional al vector 
Q(b,, db», .., by) si existe un número k tal que b,=ka,, 
b, = ka», .., b, = ka, En este caso, P =kQ. 

El concepto de la combinación lineal de los vectores es la 
generalización de la noción de proporcionalidad de los vectores. 

El vector P se llama combinación lineal de los vectores P,, 
P,, .., P,, si existen algunos números reales !,, l,, .., l, tales 
que P =!,P, + LP, + ... + 1,P,. En este caso el ¡-ésimo componente 
del vector P (i = 1, 2, ..., m) será igual a la suma de los productos 
de la multiplicación de cada ¡-ésimo componente de los vectores 
P,, P,, ..., P, por su l,, l,, ..., l, respectivo. 

Un sistema de vectores P,, P,, ..., P,_,, P, se llama linealmente 
dependiente si por lo menos uno de ellos es una combinación 
lineal de los demás vectores. 

Esta definición es equivalente a otra: 

un sistema de vectores se llama linealmente dependiente si 
existen tales números reales k,, k»>, .., k,, de los que por 
lo menos uno no sea cero, que satisfagan la igualdad 


k,P, + k,P, +... + k,P, = (0, 


En el caso contrario, el sistema de vectores se llama linealmente 
independiente. 

Si el vector P¿ es una combinación lineal de los vectores P,, 
P,, ..., P,, se dice que P, se expresa en forma lineal por medio 
del sistema de vectores (P;;, en el que ¡=1, 2, ..., n. Está claro 
que si el vector se expresa en forma lineal por medio de cierto 
subsistema del sistema dado, él se expresará también en forma 
lineal por medio del sistema para lo que es suficiente tomar los 
demás coeficientes del sistema iguales a cero. 
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Generalizando esta terminología, se dice que un sistema 


de vectores Q;,, Q,, ..., Q, se expresa en forma lineal por medio 
de un sistema de vectores P,, P,, ..., P, si cualquier vector Q;, 
i¡=1, 2, .., s es una combinación lineal de los vectores del 


sistema [Pj j=1, 2, ..., n. 
Examinemos en el espacio P“” los vectores: 
iy (1, O, O, se 0), 
i (0, 1, O, ...9 0), (1) 


i, (0, 0, 0, ..., 1). 
Estos vectores se llaman vectores — unidad (o vectores — unitarios). 


El sistema de vectores (1) es linealmente independiente, puesto que 
ki, + ki +... + Ki, =0 solamente si k,=0 para todos los 


Fig. 1 


¡=1, 2, .., m. Cualquier vector P(a,, az, ..., dm) del espacio Pl” 
se expresa en forma lineal por medio de los vectores del sistema (1), 
o sea, 


P — as1y + aL + ... + Ay - 


Se puede demostrar que cualquier sistema de vectores de un espacio 
Pp que esté compuesto por más de m vectores es linealmente 
dependiente. 

De este modo tendremos que si en una superficie salen del 
origen de las coordenadas dos vectores P, y P, que no se 
encuentran en una misma dirección, O sea que son linealmente 
independientes, cualquier tercer vector se podrá representar como 
una combinación lineal de estos vectores. 

Por analogía, si en un espacio tridimensional se dan tres 
vectores que no se encuentran en una misma superficie y que 
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salen del origen de las coordenadas, cualquier vector de este 
espacio se expresa como una combinación lineal de estos vectores. 
En la fig. 1 se ilustra el caso en el que el vector P, se representa 
por medio de una combinación de los vectores linealmente 
independientes P,, P,, P3, que tiene la notación siguiente 
Po =P; lp, +P, 
2 3 

En un espacio bidimensional, a dos vectores linealmente 
independientes P, (a,,, az) y P,(a,,, a,,) les corresponde la 
determinante 


d11 A12 
da1 22 


, diferente de cero. 


El valor absoluto de esta determinante es igual a la superficie del 
paralelógramo construido con los vectores P, y P,) (fig. 2). 
En un espacio tridimensional, tres vectores linealmente independientes 
P, (4,1, 421, 431), P¿(a,7, 422, 432) y P3(a13, 423, 433) forman 
un paralelepípedo (fig. 3). En este caso, el valor absoluto de 
la determinante será 
411 412 413 
d1 422 023 
d31 432 433 
diferente de cero e igual al volumen del paralelepípedo. 
Análogamente, si en un espacio m-dimensional se dan m 
vectores linealmente independientes 
P;(a;;, 12 j> ..., TE Ami), j= 1, Ze a IM 
O AA 


entonces, como se demuestra en los cursos de álgebra superior, 
la determinante 


0411 012 +... Olm 
431 422 ... Am 
Ami Am2 Umn 


es diferente de cero. 
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Supongamos que en un espacio m-dimensional están dados n 
vectores cualesquiera 


P;¡(41j U2j ..o, Aijo Am) A 
dl A 


Formemos con los componentes de estos vectores la matriz” 
de un orden de (m x n) 


Po 
0 p, 
Fig. 2 
P, P, P, P, 
di1 1 dy) din 
(21 422 42; 425 
o Q) 
as 
Amt  Am2 -** Ami ---  Amn 


Hablando en general, las columnas de esta matriz vistas como 
vectores m-dimensionales pueden ser linealmente dependientes. 
El máximo número de columnas linealmente independientes de 
la matriz (2) se llama rango de esta matriz**, Dicho de otra 
forma, el rango de la matriz (2) es igual al máximo número 
de vectores linealmente independientes P;,, cuyos componentes 
forman sus columnas. Cualquier sistema linealmente independiente 
máximo de vectores del espacio P” se llama base de ese espacio 
vectorial. 


*) Aquí y a continuación colocaremos los componentes 
de los vectores P;, en forma de columnas de las matrices. 
**) El número máximo de filas linealmente indepen- 


dientes de la matriz siempre es igual al número máximo de columnas 
linealmente independientes. 
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Supongamos que los vectores P,, P,, ..., P;-,, P;, Pj+1, ..., P,, 
forman la base del espacio P'"” la que llamaremos A-base. 
Entonces cualquier vector P de ese espacio puede ser representado 
en la base dada en forma de una combinación lineal de vectores 
P;, j¡=1,2, .., m. Además, esta forma de representación . tendrá 
una sola aceptación. Todas las bases de un espacio vectorial se 
componen del mismo número de vectores. 

Tomemos un vector cualquiera Q que no pertenezca a la base 
P,, P,, .., P,,. Entonces, si el vector Q no es igual a cero, 
en la combinación lineal 


Q = a: P, + a,P, O Aj-1P;- y + aP; + + ¡Pj, y A 00 E 


habrá por lo menos un coeficiente diferente de cero. 
Supongamos, por ejemplo, que a; O, entonces el vector P; 
de la A-base se puede representar de la siguiente manera 


1 0% y 1 +1 > 
P, = —OQ — —P, E TE — Po SS ZP, 1 —. — —”P. 
A A; A; A a 


Excluimos de la A-base el vector P; y añadimos a los vectores 
que quedan el vector Q. El sistema de vectores P,, P,, .... P;-4, 
O, Pj+1, ..., P,, de nuevo es base. Realmente, teniendo en cuenta 
que la descomposición de cualquier vector entre los vectores 
de la base se efectúa siempre de una sola forma, después de la 
exclusión del vector P; de la A-base, el vector Q ya no podrá 
ser representado como una combinación lineal de los demás 
vectores. Eso quiere decir que el sistema de vectores P,, P,, ..., 
Pi, Q, Pj+1, .... P, es un sistema linealmente independiente, 
oO sea, que forma una base. 

La B-base, que se ha obtenido de la A-base por medio 
de la sustitución de cierto vector por un vector no perteneciente 
a la A-base, la llamaremos base con una sustitución con respecto 
a la A-base. 

Supongamos que en un espacio m-dimensional están dados 
n+ 1 vectores Po, P,, P,, .., P,, .., P,, los que en la base 


J 
Q,, Q», ... O, tienen la descomposición 


Po, = b,Q, + b2Q) +---+ BmQm, 

P, =4,1Q, + 2210, +...+ Qm1 Om, 

P, = a12Q, + a220Q, +...+ 42 Om, 
A o (3) 


P, a 1H 01 ae 420 +h...+ Amn Qs: 
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Formemos una matriz con los componentes de la descomposición (3) 


Po P, o». P, P, 
bi 411 ... Grj ... jp 
b, da 1 e 43; ... 42 
(4) 
b; 4; a; Ain 
Dr. Und bs ds 
Admitamos ahora que los vectores Q;, Q», ..., qm forman una 


nueva base P”. Supongamos que estos vectores tienen la siguiente 
descomposición en la base Q,, Q,, .... Q,,,: 


4 = 9110, + 9210 + --- + Qm1Qm> 
d» = 91201 + 9220 + .. + 4m20m> 


. . . . .. ... ... .....o.o...o..«<o...os..o...os. ea..o. ose 


Um = dimQ 1 EN 42mQ2 ho. + ImmQ m- 


Como 4Q,, 9», -..» q, forman una base, la determinante 
411 412 -+>» Qim 
1 2 --- Um 
a / (S) 
Ami m2 --- Umm 


será diferente de cero. 

Pasemos de la base Q,, Q,, ..., O, a la base Q;, 02, ...» Qm- 
Entonces, los vectores P,, P,, P,, .., P, tendrán en la base 
Q1> 42, -..» q, Unos coeficientes de descomposición que, en general, 
serán diferentes de los coeficientes de descomposición en la base 


Q,, Q», q. O,,: 


Po = b191 + D2Q2 +... + Dindlm> 
Pp, = ar + 42 ¡2 o + Am jÚm (j = 1, Ze ad n). 


En este caso, la matriz de la descomposición tiene la siguiente 
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notación 
Po P, P; Pp. 
1 411 ay j Ar 
bz 421 a) d2n 
: (6) 
Di ar ai a; 
Í ij in 
A A A 


Como se sabe, los elementos de la matriz (6) se determinan por 
medio de los elementos de la matriz (4) y de la determinante (5) 
con las fórmulas siguientes: 


1 A ij A ( ) 
en las que 
d11 +++ Q1,i-1 by d1,i+rt «+ im 
A? =|421 --* G2ai1 Da Ga irr --» 02m 
4m1 Pia dm, i- 1 Bm dm, i+1 A Umm 
11 «++» Qrri-1 %4%rj Q1i+to +++ Qím 
A]! = 21 +++ Q2i-1 Aj 42,i+1 «+ Q2m 
4m1 cÓs Um, Y Amj Um, +1 << Umm 


=12.m  j=1,2..,m) 


A continuación, al mismo tiempo de denominar estas determinantes 
con A, A?, A?, emplearemos a veces las designaciones siguientes: 


Á = (4192 --- qm); Ay =(Q1 - qi-1P0Q;+1 ---Qm); 
A¡=(Q1 .. 4i-1PQi+1 - qm): 


Examinemos un ejemplo de determinación de la descomposición 
de un sistema dado al pasar a una base nueva. 

Supongamos que en la base P,, P,, P3, los vectores Py, P,, 
P,, P,, P,, P;, P, tienen unas descomposiciones que se determinan 
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por la matriz 


Puesto que la determinante compuesta de los componentes de los 
vectores P,, Ps, Ps 


e il 
A=lPo, P., P¿) = — 1 1 —2 = —18 
O | 1 


es diferente de cero, los vectores P¿, Ps, P¿ formarán una base. 
Encontremos la descomposición de los vectores P,, P,, P,, Ps, 
P,, Ps, P¿ en la base P,, Ps, P¿. Empleando las fórmulas (7), 


obtendremos 
Po P, P, P, 
¡ (PoP;sPó), (P,PsPó), (P>PsP¿), (P3P5P4), 
A (P¿P¿Po), (P¿P, Po), (P¿P,P,), (P,¿P3P4), 
(P,PsPo), (P¿PsP,) (P¿P,P,), (P¿PsP3), 


P, Ps Ps 
(P¿P5P4), (P¿P5P,), (P¿PsPg) 
(P¿P¿Pó), (P¿PsPo), (P¿P¿Ps) 
(P¿PsPs), (P¿PsPs), (P¿PsPs) 


Realizando las transformaciones correspondientes, obtendremos 


la matriz que representa la descomposición de estos vectores 
en la base P,, Ps, Pé: 


Po P, P, P3 P, Ps Po 


a O 
9 6 18 18 
14.11.13 0 10 
9 6 18 18 
3 124% _1l00m1 
9 6 18 18 
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Como e€s sabido, se llama producto escalar de dos vectores 
P y Q la expresión definida por la fórmula 


(P . O) = a1b, + a,b, +... + a;b; + ..+ Dos 


donde a;,, b; son los componentes de los vectores P y Q. 
Los vectores P y Q se llaman ortogonales si su producto escalar 
es igual a cero. Ya que se puede considerar a una matriz de un 
orden de mxn como un vector (m x n)-dimensional, la llama- 
remos producto escalar de las matrices 


lli1 012 +... .. ln 
ta 1 1432 e ... qa, 
Á = : ; : 
.— o... . ...,.. os. as; . . . ... 
din A Un 
y 
Di DA ate ci Di 
bz1 Daz Dan 
B= . ' 
y IR a Bic das 
Dont Dz . o. ». +... o . 0.00 Din 
a la suma algebraica de todos los elementos de la matriz 
411 Di 412 Diz ... so, A Di 
d21 Da 422 Daz... «Gan Day 
E e E la died Aids da 
Um Dis Am2 7) dá Amp ¿a 


y la designaremos (A - B). 
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$ 2, HIPERPLANO Y SEMIESPACIO 


Por la Geometría Analítica sabemos que a la ecuación lineal 
A¡X1 + Á>X> + Á3X3 E Y (8) 


en un espacio tridimensional le corresponde el plano normal 
al vector A(A;,, A), Az). 
Reduzcamos la ecuación (8) a la forma 
E EN NE A (9) 
41 42 43 
El plano que corresponde a esta ecuación, corta en los ejes 
de las coordenadas los segmentos a;,, a,, az (fig. 4). 


3 


O—origen de las coordenadas 


Fig. 4 Fig. 5 


La ecuación de una superficie plana en un espacio tridimen- 
sional puede ser, además, representada en forma vectorial 
(A. X) =h, donde A” es un vector de longitud igual a uno 
(unitario), normal al plano; X, el vector corriente que une el 
origen de las coordenadas con cualquier punto perteneciente 
al plano. Aquí (A” - X) es el producto escalar de la multiplicación 
de' los vectores A” y X que tiene un valor igual a la proyección 
h del vector X en la dirección determinada por el vector A?. 
El valor de h es igual a la distancia entre el origen de las 
coordenadas y el plano. El plano pasará por el origen de las 
coordenadas cuando h=0. 

Cualquier vector X que una el origen de las coordenadas con 
un punto del plano tendrá la misma proyección h en la dirección 
A” (fig. 5). 
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Análogamente llamaremos hiperplano, o sencillamente plano 
en el espacio m-dimensional, al conjunto de todos los puntos 


(Xi, X2, .. Xm) que satisfagan la ecuación 
AX: + ÁA>X> + ... + AL e C. (8) 
Convendremos en decir que este hiperplano es normal al vector 
A(A¡, Az, ..., Am). Diremos también que la ecuación 
A A (9”) 
a; 144) Am 


corresponde al hiperplano que corta, en los ejes de coordenadas, 
segmentos de longitud igual a a,, a,, ..., Am. 


Recta P 


Ecuación de la recta Aw=h ¿4 


Ecuación del plano 
9) AE 4%, = 33 =5 


Fig. 6 Fig. 7 


Por último, la ecuación vectorial (A”-X) =h define en el 
espacio P” un hiperplano normal al vector unitario A” que se 
encuentra a una distancia h del origen de las coordenadas. 

Una recta en un plano divide a éste en dos partes, cada una 
de las cuales se llama semiplano. En la fig. 6 se ve que una recta 
divide un plano en dos semiplanos. Añadamos que la proyección 
Oc, de los vectores X,, cuyos extremos se encuentran en uno 
de los semiplanos, son menores que h= Oc, y la de los vectores 
X,, cuyos extremos se encuentran en el otro semiplano, son 
mayores que h. 

Un plano en un espacio tridimensional también divide todo 
el espacio en dos partes, cada una de las cuales se denomina 
semiespacio. 

Por analogía, diremos que un hiperplano en un espacio 
m-dimensional divide este espacio en dos partes cada una de las 
cuales se denomina semiespacio. 
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Supongamos que un hiperplano en el espacio P se expresa 
por medio de la ecuación (A” - X) = h. Entonces, para los puntos 
M de uno de los semiespacios, las proyecciones que representan 
sus vectores X en la dirección A” serán menores que h, y para 
los puntos del otro semiespacio, mayores que h. Así diremos que 
uno de los semiespacios es el conjunto de vectores X para los 
que se cumple la desigualdad (A”-X) < h, y para los vectores del 
otro semiespacio, la (A” - X) > h. El propio hiperplano (A” - X) = h 
puede ser agregado a uno de los semiespacios. Entonces todo 
el conjunto de puntos del espacio m-dimensional se dividirá 
en dos partes: la de los puntos para los que (A”-X)<h 
y la de los puntos para los que (A-X) > h (o bien, en otro caso 
(A*-X)<h y (A”-X) > »). 

Ejemplo 1. Dado el semiespacio —5x, + 4x, — 3x3 < 5 deter- 
minar si le pertenece el punto (0, 0, 0). 

Para obtener la respuesta es suficiente colocar los valores 
x1=0, x=0, x3=0 en la desigualdad. Tendremos —5-0+ 
+4-0-3-0<5, de lo que se deduce que el punto (0, 0, 0) 
pertenece realmente al semiespacio —5x, + 4x, — 3x3 < 5. 

El plano —5x,+4x,-—3x3=5 es perpendicular al vector 
A(—S, 4, —3) (fig. 7). 

Ejemplo 2. Determinar si el punto de un espacio 9-dimensional 
x1=0, x2=4 x3=3, x4=-—7, x5=0, x=0, x7=9, x3=1, 
xy =0 pertenece al semiespacio 4x, + 5x, — 7X3 + x5 — 2x6 + 
+ 12x, — 3x9 < 11. Colocando las coordenadas de este punto en la 
desigualdad, obtendremos 4-04+4-5-3.7-7-0+0-1-2-0+ 
+9.12+0-1-3-0= 107 > 11. En consecuencia este punto no 
pertenece al semiespacio dado. 


$ 3. POLIEDROS CONVEXOS 


Se llama cuerpo convexo a un cuerpo que contenga junto con 
dos cualesquiera de sus puntos, todo el segmento rectilineo que los 
une. 

Ejemplos de cuerpos convexos son un círculo, una esfera, 
un cubo, un ángulo formado por dos líneas que salen de un 
punto (fig. 8), etc. 

Supongamos que los puntos x e y son comunes para los 
cuerpos convexos A y B (fig. 9). Entonces x e y pertenecen al 
cuerpo ÁA y por eso el segmento que une los puntos x e y 
también pertenece a A. Por analogía, este mismo segmento 


Zal 


pertenece también al cuerpo B. Por consiguiente él también 
pertenece a la parte común de los cuerpos A y B. Esto significa 
que la parte común o la intersección de los cuerpos convexos es 
un cuerpo convexo. 

En una superficie plana tomemos un polígono que se encuentre 
siempre a un lado de cada recta que forma este polígono. Este 
polígono, como se ve en la fig. 10, a, es convexo. Realmente, 
dos puntos cualesquiera x e y de este polígono le pertenecen 
junto con el segmento que les une. Al contrario, el polígono 
representado en la fig. 10, b no se encuentra por entero en un lado 
de cada una de las rectas que forman este polígono. Un polígono 
tal no es convexo. 


Fig. 8 Fig. 9 


Tomemos un punto cualquiera M que no pertenezca a un 
polígono convexo (fig. 11, a). Siempre se puede trazar una recta PQ 
tal que el punto M y el polígono queden a diferentes lados 
de ella. Se puede construir para un polígono convexo un conjunto 
de rectas tales que cada una de ellas tenga por lo menos un punto 
común con el polígono y que todo el polígono quede a un lado 
de cada una de ellas. Estas rectas se llaman soportes. 

En el jemplo de la fig. 11,b las rectas AD, CB, DB y FN 
son soportes. La recta soporte puede tener una parte común 
con el polígono convexo. Esta puede ser un punto o un segmento. 

En un espacio tridimensional, un cuerpo limitado por planos 
que se encuentra por entero a un lado de cada plano que 
contiene una de sus caras, es convexo y se llama poliedro convexo. 
Ejemplos de tales cuerpos son un diamante poliédrico, un prisma, 
etc. El plano que tiene por lo menos un punto común con el 
poliedro y le deja a todo él a un lado de si se llama plano 
soporte de un poliedro convexo. El plano soporte puede tener con 
el poliedro una parte común que puede ser un punto (llamado 
vértice del poliedro), un segmento (llamado arista), o un polígono 
(llamado cara). Está claro que por cada vértice o arista del 
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poliedro se puede trazar una cantidad infinita de planos soportes, 
mientras que por cada cara puede pasar sólo uno. 
Anteriormente nos convencimos de que la parte común 
de varios cuerpos convexos es un cuerpo convexo. Por eso, 
la parte común de varios poliedros convexos es también un cuerpo 
convexo. Ya que el plano es un cuerpo convexo, la intersección 
de un poliedro con un plano también es un cuerpo convexo que 


Fig. 11 


puede ser un punto, un segmento o un polígono convexo. 
Las propiedades de los cuerpos convexos de los espacios multi- 
dimensionales se pueden examinar análogamente a las propiedades 
de los cuerpos convexos del espacio tridimensional. Algunas 
de estas propiedades se estudiarán en los $$ 4 y $. 


$ 4. SISTEMA DE DESIGUALDADES LINEALES 


Supongamos que en un espacio bidimensional se dan n 
desigualdades de la forma 

A¡1X1 + Ai2X < b; (i = 1, 2, es n) *) (10) 

*) Cualquier desigualdad de la forma a;,x, + 4¡2X> > D), 


después de multiplicar sus dos miembros por —1 se convierte en la 
forma (10). 
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Cada una de estas desigualdades determina uno de los dos 
semiplanos que tienen por línea límite la recta a;¡x, + ay,Xx, = b;. 
La recta-límite a;,x, + ax, = b; es normal al vector A, (a;,, a;»). 


xo 


Ecuación de recta 
2x+3x56 


Ecuación de recta —> 
+A FE 


a) b) 


Fig. 12 


Llamaremos solución del sistema dado a cualquier par de 
números (x,, x2) que satisfaga todas las desigualdades del sistema 
(10). En otras palabras, cualquier punto del plano (x,, x,) cuyas 
coordenadas satisfagan el sistema (10), es una solución. 
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Veamos algunos ejemplos: 
1. La desigualdad 


X1 


3 


X2 


A Ó 2x1 +3x,<6 


determina un semiplano (fig. 12, a). A esta desigualdad le satisface 
cualquier punto que se encuentre en la parte rayada del plano. 
La recta-límite se expresa con la ecuación 2x,+3x,=6 y es 
normal al vector A (2, 3). 

2. Dos desigualdades 


2x1 + 3x2<6 y —=Xx1+x,<2 


determinan una parte del plano, como se muestra en la fig. 12, b. 
3. A las tres desigualdades 


les satisface el conjunto de puntos del plano que forman el 
triángulo de las soluciones AFB (fig. 12, c). 
4. A las cuatro desigualdades 


2x1 + 3x, <6, (1) 
—X1 | X2 < Ze (2) 
—xXx, — 3x, < 3, (3) 

3 
X1 5 (4) 


les corresponde el conjunto de puntos del plano que forma 
el polígono de las soluciones ABCD (fig. 12, d). 
5. A las siete desiguladades 


2x1 +3x <6, (1) 
—=X1+xXx <2, (2) 
—x1 —3x, <3, (3) 
DN <= 3, (4) 
201 < 3, (5) 
—3x, +%7x, < 21, (6) 
xi — 3x, < 3 (7) 
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les corresponde el mismo conjunto de puntos que en el ejemplo 4. 
Las desigualdades (5), (6) y (7) pueden ser excluidas sin alteración 
del conjunto de soluciones. En este caso, las desigualdades (5) y (6) 
definen las rectas-límite que no tienen puntos comunes con el 
polígono ABCD. La recta (7) tiene un punto común con 
el polígono y es soporte (fig. 12, e). 

6. El sistema de desigualdades 


2x1 + 3x2 <6, (1) 
—=X1+X2 <2, (2) 
—Xx1 — 3x2 $3, (3) 
2x1 0 (4) 
—3x, — 2x, < — 12, (S) 


no tiene ni una solución. Desde el punto de vista de la Geometría, 
esto quiere decir que no existe ni un punto cuyas coordenadas 
satisfagan todas las desigualdades (fig. 12, f). 

El estudio de estos ejemplos nos lleva a las siguientes 
conclusiones: 

1. Un sistema de desigualdades con dos variables puede ser 
compatible. Entonces existe por lo menos un punto del plano que 
pertenece a todos los semiplanos determinados por el sistema dado 
de desigualdades. El conjunto de todos los puntos tales puede ser 
un semiplano, un polígono limitado o ilimitado*, una recta 
o su segmento y, por fin, un punto. El conjunto de puntos que 
satisfacen el sistema de desigualdades es un cuerpo convexo. 

2. Un sistema de desigualdades puede ser incompatible. En este 
caso, no existe ni un punto del plano que satisfaga simultáneamente 
todas las desigualdades del sistema. 

En un espacio tridimensional, un sistema de n desigualdades 
puede escribirse, sin limitarle la generalidad, en la forma 


411 X1 + 42X2 + Mi3X3<b; (i=1, 2, ..., n). (11) 
Como ya sabemos, cada una de las desigualdades (11) define 
un semiespacio con el plano límite 


41X1 + 0¡2X2 + M¡3X3 = Dj. 


*) Pueden haber sistemas de desigualdades a los que 
correspondan figuras (poliedros) no limitados por sus aristas y caras en 
alguna de las direcciones. Estos sistemas de desigualdades también 
responden a esta propiedad. 
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El conjunto de desigualdades (11) puede ser compatible. En este 
caso, existe cierto conjunto de puntos del espacio tridimensional 
que satisfacen el sistema de desigualdades. Este conjunto es un 
conjunto convexo y puede representarse en forma de un semiespacio, 
un poliedro, un plano, un polígono, una recta o, por último, 
un punto. 

En el caso de la compatibilidad, entre las desigualdades del 
sistema puede haber desigualdades “excesivas”, que al ser eliminadas 
no modifican el conjunto de las soluciones del sistema de 
desigualdades dado. Estas desigualdades “excesivas” pueden ser 
de dos clases. Las de la primera clase son aquellas cuyos planos- 
límite no tienen intersección con el conjunto de todas las soluciones 
del sistema. Las de la segunda clase son las desigualdades cuyos 
planos-límite son soportes para el conjunto de soluciones. 

Puede resultar que en un espacio de tres dimensiones no haya 
ningún punto que satisfaga simultáneamente todas las desigualdades. 
En este caso, el sistema de desigualdades se llama incompatible. 

Supongamos que en un espacio m-dimensional se da el sistema 
de desigualdades 


AX + Op9X2 +... + QimXm S bj, (i=1, 2, ..., n) (12) 


Por analogía con el espacio tridimensional diremos que cada 
una de las desigualdades (12) define en el espacio m-dimensional 
un semiespacio con el plano-límite 


A¡1X1 + A2X2 +... + UimXm = b;. 


Si en el espacio m-dimensional existe por lo menos un punto 
M (Xi) X2,. ... Xm) que satisface simultáneamente todas las 
desigualdades, el sistema (12) se llama compatible. Al conjunto 
de todos los puntos tales lo llamaremos “el poliedro de las 
soluciones”. 

Supongamos que en un espacio m-dimensional se dan dos 


puntos M'(XL, X2> .«.» Xm) Y M'"(Xi, X2, ..., Xm) El conjunto 
de puntos M(X;¡, X, +...» Xm) cuyas coordenadas satisfacen las 
condiciones 


X1 =X1 +1 (xi — x1), 


x2 + 1 (x2 — X2), 


e 
ES 
l 


Xm FE Qin — Xm) 


e 
S 
l 
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al variar el parámetro t de O a 1, se llama segmento que une los 
puntos M' y M”. El “poliedro de las soluciones”, siendo la 
intersección de los semiespacios, es un conjunto convexo. Eso quiere 
decir que junto con los puntos M' y M” también le pertenecen 
todos los puntos del segmento que los une. 

Las desigualdades que pueden ser eliminadas del sistema (12), 
sin modificar el conjunto de sus soluciones, se llaman dependientes 
o “excesivas”. Al eliminar consecutivamente del sistema dado una 
desigualdad de esta clase tras otra, obtendremos un subsistema 
de desigualdades con un conjunto de soluciones que coincidirá 
con el conjunto de soluciones del sistema inicial. 

La eliminación de las desigualdades “excesivas” es un proceso 
muy complicado y trabajoso. Una de las particularidades de los 
métodos de programación lineal consiste en que para calcular 
el valor mínimo (máximo) de la función lineal en el poliedro 
no es imprescindible eliminar las desigualdades “excesivas”. 
Si al sistema de desigualdades (12) no le satisface ningún punto 
del espacio m-dimensional, este sistema se llama incompatible. 


$ 5. EL VALOR MÍNIMO Y MÁXIMO 
DE LA FORMA LINEAL EN EL POLIEDRO 


Examinemos un sistema compatible de ecuaciones lineales con 
dos variables. Consideremos que ya hemos eliminado todas las 
desigualdades “excesivas” de la 1-a y de la 2-a clase, y así hemos 
seleccionado el polígono de las soluciones en su “forma pura” 
(fig. 13) Supongamos que, además, se da la función lineal 
de dos variables 


f = C1%X1 + C2X2- 


Encontremos en el conjunto de puntos (x,, x,) del polígono 
de las soluciones tales puntos que lleven la función lineal 
f =c,x, + c,x, al valor mínimo y máximo. Examinemos el conjunto 
de todos los puntos (x,, x,) del plano en cada uno de los 
cuales la función f =C,x, +C,x, toma un valor fijo (f =f,). 
El conjunto de tales puntos es la recta c,x, + c,x, = f,. Esta recta, 
como se ha señalado en el parágrafo anterior, es normal al vector 
C(c,, c,) que sale del origen de las coordenadas. Tracemos una 
recta F (fig. 13) normal al vector C y desplacémosla paralelamente 
a sí misma en la dirección positiva del vector C. Supongamos que 
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en su desplazamiento, la recta F se encuentra por primera vez 
con el polígono en el vértice A. En esta posición F' la recta F 
se hace soporte. Al continuar el desplazamiento en esa misma 
dirección, la recta F pasará por el vértice B y se hará también 
recta-soporte. Puesto que el sentido del vector C(c,, c,) es la 
dirección del máximo incremento de la función lineal f =c,x, + 
+ cx, entonces entre todos los valores que toma la función 
lineal f en el poligono de las soluciones, esta función tomará 
en la recta-soporte F' su valor mínimo y en la recta-soporte F”, 
su valor máximo. 


Fig. 13 Fig. 14 


Así que los valores mínimo y máximo de la función lineal 
f =c,x¡ +C€,x, en el poligono de las soluciones se alcanzarán 
en los puntos de intersección de este polígono con las rectas- 
soportes normales al vector C(c,, c,). La intersección de la recta- 
soporte con el polígono de las soluciones puede ser un solo 
punto* (un vértice del polígono) o un conjunto innumerable 
de puntos (en este caso el conjunto es un lado del polígono). 

En la fig. 14 se presenta el caso en el que la función lineal f 
alcanza el valor mínimo en cada uno de los puntos del segmento 
A¡,, Az, mientras que el valor máximo, en los puntos del polígono 
que están infinitamente alejados. 


*) Puede resultar que este punto se encuentre en el 
infinito. 


3) 


Por analogía, una función lineal de tres variables f =c,x, + 
+ C2X + C3x3 toma un valor constante en un plano normal 
al vector C(c,, c», C3). El sentido de C es la dirección del máximo 
incremento de la función f (fig. 15). Los valores máximo y mínimo 
de esta función en el poliedro de las soluciones también se alcanzan 
en los puntos de intersección de este poliedro con los planos-soportes 
normales al vector C(c,, c,, ca); aquí la función f alcanza en uno de 
los planos-soportes el valor mínimo y en otro, el máximo. La inter- 
sección de un poliedro con un plano-soporte puede ser un punto 


X3 


Fig. 15 


Poliedro de las 


Soluciones es 
O 


E 
Plano 
soporte 
A Plano 
sopor te 


Fig. 16 


(un vértice del poliedro) o un conjunto innumerable de puntos 
(este conjunto será un arista o un lado del poliedro). 

Por ejemplo, en la fig. 16 aparece el caso en el que f alcanza 
su valor mínimo en cada uno de los puntos del lado ABOP 
y su valor máximo, en el punto D. 

3* 
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La generalización del concepto de las funciones lineales de dos 
o tres variables es la función f =C,x,+CC>,Xx3+..+C,x, de n 
variables reales x,, X», .., Xp, llamada forma lineal. Aquí c,, 
Ca, «.» C, SON números reales. Fijando valores de la forma lineal, 
o sea f=f13 f=f23 .; f=f, se definen en un espacio 
n-dimensional los hiperplanos 


f, = C1X1 + C2X +... + CrX y 
f, =C,X1, + C2X2 + ... + CyXp 


f¿ =C1X1 + C2X2 + ... + CpXp 


que son normales al vector C(c,, Cz, .., Cn). Designemos con 
f mn el valor mínimo de la forma 


f = C1X1 + C2X +... + CrXn 


en el poliedro de las soluciones, y con fix, el valor máximo. 
Puesto que el vector C determina la dirección del incremento 
máximo de la forma lineal f, entonces cuando las funciones 
FPS Fmin Y > foma, los hiperplanos correspondientes f' =C,x, + 
+ C9X2 +... + C¡X, no tienen intersección con el poliedro de las 
soluciones. Por otra parte, cada hiperplano f” =C,x, + C,X, + 
+... + CpX,, para el que f min € f” € f máx, tiene puntos comunes 
con el poliedro de las soluciones. El conjunto de puntos M (x;, 
X2» ...» Xp) en los que la forma lineal alcanza el valor mínimo 
es la intersección del poliedro de las soluciones con el hiperplano 
SOPorte C¡X¡ + C2X2 +... + CrXn = f min, normal al vector C(c;,, 
Cz» «. Cy). Análogamente, f alcanza el valor máximo en los puntos 
de intersección del poliedro con el hiperplano-soporte c,x;, + 
+ C7X2 +... E CaXn = fm, que es también normal al vector C. 

La intersección del poliedro de las soluciones con un 
hiperplano-soporte será un vértice, una arista o un “lado” del 
poliedro. 

Así que el valor óptimo de la forma lineal en el poliedro 
de las soluciones se alcanza en los puntos entre los cuales siempre 
se encuentra aunque sea uno de los vértices. Por eso para 
calcular la solución óptima es suficiente encontrar el vértice del 
poliedro en el que la forma lineal alcanza el valor mínimo 
(máximo). 
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$ 6. REDUCCIÓN DE DESIGUALDADES A IGUALDADES 
AL RESOLVER PROBLEMAS DE PROGRAMACION LINEAL 


Se da un sistema de m desigualdades lineales con n variables: 


111X1 + 012X2 + -.. + 01¡X, < 
1421%X1 + 4A22X2 +. + A2nXn < 


O 5 (12) 


. . . . . . .oxo.no..ox2.02..kn.02x2.P.02o2.2.oxÁ2b.2.b2xXoxÁnbnbxXx—26%xxX—b_hb . .<.. .. .. 


AmiX1 + Am2X2 +... + AmnXp < ¿ON 


que determina en un espacio n-dimensional un poliedro de las 
soluciones. Simultáneamente a este sistema de desigualdades veamos 
un sistema de m ecuaciones algebraicas lineales con n + m variables: 


411X1 + 012X2 +... + 04pX, + M1 =D, 
A21%X1 + 1d22X2 +... + A2nXn + y = b,, 
tobas ana (12) 


. . . . .a.... . .oa.oo. .ob..2o.0202...20.o2.o2...2P2...—.—_. ..o..o.2.22.02.0202... 


Am1iX1 + Am2X2 +. + AmpXy + Mm = Bin 


Demostremos que a cualquier solución xi, x23, .... x; del sistema 
de desigualdades (12) le corresponde una determinada solución 
Xi» Xd, «0. Xp3 Hi, 43, »., Hp del sistema de ecuaciones algebraicas 
lineales (12”), teniendo en cuenta que las variables complementarias 
(variables de holgura) satisfacen la condición de que uz > 0, 1% > 0, 
..., Hp > 0. Realmente, si xi, X2, ..., Xx, es la solución del sistema 
(12), tienen lugar las desigualdades: 


ay (X1 + A412X2 a (11 yA 
471X1 + 079X3 +... + 02,X 


. . .. o. +. .—...o.2. . .2...2.2.020.ox.>£2>._2.o02.2x. . . ss.» 


. . .. .. .. .. . .. ...o. .o..—...ea<.oxÁX£n20 2020... o. o 


RO 1 A E E 1 a y 
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Designemos 


Le] 


a] o o 
pi =b, —(84,,Xx1 + 412X3 +... + 0,,Xp), 


(0) O 10) 10) 
3 = b — (471X] + 492X3 +... + 02yXp), 


A A dl a A o A A (12) 


- o o o 
Um = Br aa (Am 1 X1 LE UAm2*X2 Ed a "o es 


Entonces, en primer lugar, 1í > 0, 13 >0, ..., 4, > 0; en segundo 
lugar, como se deduce de (12”) el sistema de números xi, X2, ..., 
Xp3 Hi, H3, ... Hp es la solución del sistema de ecuaciones 
algebraicas lineales (12). 

Demostremos, a la inversa, que a cualquier solución xi, X2, ..., 
Xp3 Hi, 45, ..., Hp del sistema de ecuaciones (12), que satisfaga 
la condición 13 > 0, 13 > 0, ..., uy, > 0, le corresponde una solución 
determinada del sistema de desigualdades (12). En efecto, como 
el sistema de números xi, X2, .... Xp3 HI, 43, «.., Hp €s la solución 
del sistema (12), podemos escribir: 


lo) O o) o) 
411%X1 + 412X23 + 27s + A1nXn + Hg — b,, 


._. . . .. . .....o.o..s..o..o o. o ..—..osoeose 


. . .. .. ..........o..o..o..o..o.o.o.o..os.o.o.o.o.o se eo ss 


0) 9) (0) 
AmiX1 e Am2X2 EN AmnXn ES Mn > Dis 


De acuerdo a las suposiciones hechas los números Hi, H3, +...» Hp 
son no negativos, entonces se cumplen las desigualdades: 


(0) [0] (0) 
4114X1 + 019X2 +... +01,X, <bi, 


0) 0) O 
AmiX1 + Am2X2 +... + AmnXn < Dm, 


o sea, el sistema de números xi, X2, ..., Xx, es la solución del 
sistema de desigualdades (12). 

De este modo se ha establecido la existencia de una relación 
recíproca entre el conjunto de todas las soluciones xi, X3, ..., Xy 
del sistema (12) y el conjunto de las soluciones xi, X2, ..., Xp; 
Hi, H2, ..., HU del sistema (12') en los que se encuentran los mismos 
valores. Observemos que las variables complementarias Hi, 42, ..., 


Hs, son no negativas. O sea, el problema de la resolución de un 
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sistema de desigualdades lineales (12) se reduce a la resolución 
de un sistema correspondiente de ecuaciones lineales (12”. 

En los problemas de programación lineal mos interesarán 
soluciones de sistemas de desigualdades que satisfagan la condición 
x1>0,x2>0, ..., x, > 0. Estas soluciones se llaman no negativas. 
Por eso, como la resolución de un sistema de desigualdades se 
reduce a la resolución de un sistema correspondiente de ecuaciones 
lineales algebraicas, una de las tareas más importantes de la 
programación lineal es el cálculo de las soluciones no negativas 
de un sistema de ecuaciones lineales. 

Ejemplo. Calcular la solución no negativa del sistema de 
desigualdades: 


Introduciendo las variables complementarias u, > 0, 4, > 0, uz > 0, 
obtendremos el sistema de ecuaciones: 


2x1 + 3x2 + Hi = 6, 
+ +m o =2 
—X1 — 3x, + Hz = 3. 


A cualquier solución no negativa de este sistema de ecuaciones 
le corresponde una determinada solución no negativa del sistema 
inicial de desigualdades y viceversa. Por ejemplo, a la solución 
no negativa xy =1l, x2=1l, 4 =1, 1,=2, 143 =7 del sistema 
de ecuaciones le corresponde la solución no negativa x, = 1, x, = 1 
del sistema de desigualdades inicial. 

Observemos que si el sistema de desigualdades es incompatible 
en la región de las soluciones no negativas, el sistema correspondiente 
de ecuaciones lineales no tiene ni una solución no negativa. 

Es sabido que el valor mínimo (máximo) de la forma lineal, 
en un poliedro determinado por un sistema de desigualdades, 
se alcanza en un determinado vértice del poliedro de las soluciones. 
Se puede demostrar que a cada vértice del poliedro de las soluciones 
le satisface una solución no negativa del sistema correspondiente 
de ecuaciones algebraicas lineales en el que por lo menos una 
variable complementaria es igual a cero. Las variables comple- 
mentarias 1, > 0, 1, > 0, .... Um > O que se introducen en el sistema 
de desigualdades se pueden explicar en un sentido geométrico. 

Examinemos el poliedro de las soluciones del sistema de 
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desigualdades: 
411%Xi + A12X2 +. + Ar nXn < bi; 


Am1iX1 AN Am2%X2 os AmnXn < Bm 


Supongamos que existe una de las soluciones no negativas xí, 
X2» «.., Xy del sistema de desigualdades. A esta solución le corresponde 
el punto M(xi, X2, +... X;5) que pertenece al poliedro de las 
soluciones. 

Supongamos que con h; se denomina la distancia entre el punto 
M y el i-ésimo hiperplano. Los valores de las variables comple- 
mentarias 


o 


pj =b, — (4, 1x1 + 412X2 +... + 417Xp), 
Min = Py — (Am1X] + Om2X2 + «<< + Gn Xn) 
serán proporcionales a las distancias entre el punto M(x1, X3, ..., 
Xx.) y los hiperplanos-límites: 
411X1 + A12X2 +... + 0A19Xn = bs, 


AmiX1 + Am2X2 +2 + AmnXn = Dip - 


Realmente, según la definición, la distancia entre el punto M 
y el ¡-ésimo hiperplano es igual a 


his b, — (411Xi + Oj2X2 +... + GjnXp) 
Vaz + 4% +... + Aía 
(= 1:02: 305 10), 
de lo que se deduce que: 
y =h, Y az, E 


Mm = Am Y Ar + din2 +. + Gimp: 
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CAPÍTULO 1 


RESOLUCIÓN DEL PROBLEMA GENERAL 
DE LA PROGRAMACIÓN LINEAL 


La programación lineal es una disciplina dentro de las 
matemáticas en la que se estudian los métodos de resolución del 
siguiente problema: se da un sistema de m ecuaciones algebraicas 
lineales que contienen n variables: 

0411X1 + 0O129X2 +. + ad1¡X; +... + 01pX, = b,, 


421X1 + 4A22X2 +..+ d7¡X; +... + A2yXp = b,, 


A A E A IR (13) 
AX + OjX2 4h .. + GjjXj +... + AjynXp =D; 
Ami X1 + Ama2X2 bo + AmjX; +... + AmnXn = bn 
y la forma lineal f de estas mismas variables 
f =C0,/X1 + C2X2 +... + CjX¡ Hb ... + CpX, (14) 


Hay que encontrar entre las soluciones del sistema (13) una 
solución no negativa en la que la forma lineal f tome su valor 
mínimo *. A esta solución la llamaremos la solución óptima 
de este problema. 

En este capítulo se estudia un método de cálculo de una 
de las soluciones no negativas del sistema (13) y uno de los 
algoritmos de resolución del problema general de la programación 
lineal. 


$ 7. TRANSFORMACIONES IDÉNTICAS DE 
UN SISTEMA DE ECUACIONES ALGEBRAICAS LINEALES 


Se da un sistema (13) que contiene r ecuaciones linealmente 
independientes y resuelto con relación a r variables. Sin limitar 
la generalidad, se puede considerar que el sistema está resuelto 


*) El problema del cálculo de una solución no negativa 
en la que f obtiene su valor máximo se formula en forma análoga. 
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respecto a las primeras r variables de las primeras r ecuaciones: 
X] = by — (Air+1Xp+1 + OA1r+2Xp+2 a A1nXn), 


X2 = b2— (A2r+1iXr+1 + M2r+2Xp+42 5. + O2pXp), 


E De (15) 


. . . . .. .o...s..o...s.s.....os. .oo. ...o.o.o..os.o<ss.<s ose. 


Xp = b, ES (Ar + ra E OA 2 Xp 42 E o.» E Ay X a). 


En una anotación abreviada el sistema de ecuaciones (15) toma 
la forma 


X¡ = b; eS de 00%) | Ze da, YO (15) 
j 


Supongamos que los miembros libres b; en las ecuaciones (15) son 
no negativos. 

Cada una de las ecuaciones del sistema (15) se puede conceptuar 
como una proyección de la ecuación vectorial 


Y x¡P; = Po — y sp) (16) 
i=1 j=r+1 

en los vectores P,, P,, .., P,, donde P, =P, (1, 0, ..., 0); P, = 
= P, (0, 1, ..., 0) ..; P, =P, (0, 0, ..., 1). Los vectores P,, P,, ..., P, 
forman una base en un espacio r-dimensional. La matriz de 
la descomposición de los vectores P,, P,, P,, ..., P,, ..., P, en la base 
P,, P,, .., P, se representa en la forma siguiente 


Pi Es se Bposia Bes o ya e E P: 
1 ' E br Arr+ro. d1j + 01, 
pb E ba Oar+1 «Gay. Gay 
(17) 
1 . Pi Oir+1 ST Ain 
A AN An 


Admitamos que el conjunto de vectores P,, P,, ..., P,, P,,;,, ..., P, 
que entran en la ecuación (16) forma varias bases de un r-ésimo 
orden que se diferencian una de otra aunque sea en un vector. 
Llamaremos a cada una de estas bases base propia. 
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Sean A,, A, ..., Ay las bases propias de un sistema de vectores 
P,, Pz, .... Pj, ... P,; Prri) «o. Pj, .... P,. Al pasar de una base 
a otra, todos los coeficientes de las ecuaciones (15) se transformarán 
según las fórmulas (7) del capítulo anterior. Entonces el sistema 
(15) se convertirá en sistemas equivalentes a sí mismo o sea, 
en sistemas con las mismas soluciones. 

Llamaremos la base A positiva con relación al vector cero Po, 


si P¿ se representa en esta base en la forma 


r 
Po= B;P,, B; > 0. 

¡Z 
Como se deduce de las ecuaciones (15), la base A que contiene 
los vectores P,, P,, ..., P,, correspondientes a las variables resueltas, 
es positiva, puesto que según la suposición b;>0, Las variables 
(X1, X2, . X,), respecto a las cuales el sistema (15) está resuelto, 
se llaman básicas; las restantes, no básicas. La solución de un 
sistema que se obtiene al igualar a cero las variables no básicas 
se llama solución básica. 

En virtud del teorema sobre la descomposición de un vector 
de una sola forma en cada base, la solución básica se determina 
de una sola forma en cualquier base propia. 

A continuación vamos a estudiar los sistemas de ecuaciones 
algebraicas lineales solamente en bases propias positivas. Por eso 
para reducir la escritura omitiremos las palabras “propia” 
y “positiva”. 

Examinemos el sistema (15). Definiremos las reglas para pasar 
de una base a otra. Para ello emplearemos la matriz (17). 

I. Encontramos una de las columnas ¡(r +1 <j¡<m) en la cual 
entre los coeficientes a;;, 47;, ..., Ajj, »«.» Ap; (1 <i< r) de la variable 
no resuelta x;, hay coeficientes positivos*? (en el $ 8 se examina 
el caso en el que no existe una columna tal). 

Supongamos que ésa sea la columna j;. 


Il. Determinamos mín (e por todos los ¡ para los cuales 
Gi, 

son positivos. 

Supongamos que una de las relaciones mínimas sea aquella 

en la que i =1,. Llamaremos al elemento ai,j, resolvente (o pivote) 

con relación al sistema (15). 


di 


*) Se toman en cuenta los signos de los coeficientes 
que se encuentran dentro de los paréntesis en la notación de las 
ecuaciones en forma del sistema (15). 
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III. Resolvemos la ¡,-ésima ecuación con relación a la variable 
Xx;, y la sustituimos por su expresión en las demás ecuaciones del 
sistema (15). Entonces el sistema (15) se transformará en el sistema 
(18) que le es equivalente: 


bi, Xi, N 
xy = e was xp 
Gu dd E Gi, 


(1 = 1, Ze ..., ls —_ 1, li + 1, E 1), 


b, A»; a; 
Es 5% PU 2 E A ae YY 
s=(» pS a e A 
hi 111 NN 


¡IF 
G=r+1, ..., Mm, 
y el grupo de vectores P,, P,, ..., P, _,, P; , P; ,¡, ..., P, se convertirá 
en la nueva base, 4,. Realmente, del 18) se deduce que estos 


vectores forman una base, puesto que la determinante que les 
corresponde es diferente de cero 


Pi Bo Pio Pl ar ssp E 
1 My, 
1 A», 
1 = 4, >0 
AA 
1 
a; 1 


Esta base es positiva porque los miembros libres nuevos son no 
negativos. Efectivamente, los miembros libres tienen la forma 


b; , 
> y el miembro 
Tr 


ij, 
libre será no negativo; si a; < 0, entonces el miembro libre de nuevo 
será no negativo. 
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La transformación del sistema (15) en el sistema (18) determinada 
por las reglas 1, II, III, se llama transformación identica o simplex. 

De la definición de transformaciones idénticas se deduce que 
un sistema de ecuaciones algebraicas lineales del tipo (15) puede 
ser sometido cada vez a estas transformaciones si en las partes 
derechas hay coeficientes positivos de las variables. La sucesión 
de las transformaciones idénticas de un sistema de ecuaciones 
algebraicas lineales puede ser interpretada geométricamente como 
una descomposición sucesiva de la ecuación vectorial correspondiente 
a este sistema en las diferentes bases propias positivas que aparecen 
como resultado de una sola sustitución cada vez (o sea en cada 
transformación). 

Teniendo en cuenta que las transformaciones idénticas se emplean 
para la resolución del problema fundamental de la programación 
lineal, las estudiaremos con la siguiente condición complementaria: 

Supongamos que la última r-ésima ecuación del sistema (15) 
permite una sucesión infinita de transformaciones idénticas de este 
sistema *. Eso quiere decir que cualquiera que sea el número 
de transformaciones idénticas que hagamos, en la r-ésima ecuación 
siempre habrá coeficientes positivos de las variables en su parte 
derecha, y en cada transformación, el elemento resolvente no 
pertenecerá a esta ecuación. 

A la sucesión de transformaciones idénticas le corresponde 
la sucesión A, —> A, > 43... de los pasos de una base a otra. 
En el caso de que r y n sean finitos, existe un número finito 
de diferentes bases; entonces en la sucesión 4, > 4, > Ay >... debe 
haber regreso a bases que ya aparecieron. Diremos que en este 
caso se forma un ciclo cerrado de los pasos de base a base o que 
tiene lugar la “periodicidad”. 

Si en cada transformación idéntica la menor relación de 


Ai, 


Aj, donde j, es el número de la columna que contiene el elemento 
pivote a,, es diferente de cero, entonces el valor del miembro 


b. 
mín q , miembros libres hb; a los correspondientes coeficientes 


libre en la r-ésima ecuación del sistema (18) a cada paso sólo 
disminuye, lo que se percibe directamente de la expresión de este 
miembro libre. La elección de la base determina de una sola forma 
el valor de los miembros libres del sistema de ecuaciones, por eso 


e 


*) Los razonamientos que van a continuación son 
justos para cualquier ¡-ésima ecuación del sistema (15) (1 < ¡< r). 
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a diferentes bases les corresponden diferentes conjuntos de miembros 
libres y viceversa. Por esta razón en el caso dado es imposible 
el regreso a una de las bases que se obtuvieron anteriormente. 
En consecuencia, no puede haber una sucesión infinita de 
transformaciones. Esto quiere decir que después de un número 
finito de transformaciones, o bien todos los coeficientes de las 
incógnitas de la parte derecha de la r-ésima ecuación se hacen 
no positivos, a,¡<0 (r+1<j<nm), o bien el elemento pivote 
en cierta transformación resulta perteneciente a esta ecuación. 


S1 comenzando por cierta transformación, el elemento pivote 
perteneciese a ecuaciones con el miembro libre b;, nulo, el valor 
del miembro libre en la r-ésima ecuación permanecería invariable. 
En efecto, en la última ecuación del sistema (18) se ve que en este 


caso es igual a cero el sustraendo b, — a,;, a de la expresión 
1 

del miembro libre. En estas condiciones surge la posibilidad de 

la periodicidad. 

Llamaremos hiperplano básico al hiperplano que pasa por 
lo menos por dos vectores básicos. 

La existencia de miembros libres nulos en las ecuaciones del 
sistema que se transforma, significa que el vector P, pertenece 
a uno o a varios hiperplanos básicos. Siguiendo la terminología 
aceptada en la literatura, diremos que en el último caso tiene lugar 
una situación de degeneración. 

Asi que la condición necesaria para que ocurra la periodicidad 
es la existencia de la degeneración. 

Veamos, en calidad de ejemplo, el sistema de ecuaciones: 


X1 — Z Ss (2x4 E Xs), 
x2 =3 — (3x4 — Xs), 
X3 = 2 q (Xa + 2x5). 
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La ecuación vectorial correspondiente tiene la forma x,P, + 
+ x2P, + x3P3 = Po — (x4Pa4 + x5P5), en la que los componentes 
de los vectores determinan la matriz 


P, Pz P3 Po P, Ps 
1.0.0. 2. 2. 1 
0. 1.0. 3 3 -1 
0.0. 1. 2 1. 2 


Los vectores P,, P,, P3, P,, Ps (fig. 17, a) forman varias bases, 
por ejemplo, las bases P,, P,, P3, y P,, P3, P,. A la base P,, 
P,, P, le corresponde la solución básica x, =2; x, =3; x3=2; 
X4 =0; xs =0. En este caso el vector P¿ es una combinación 
lineal con «coeficientes positivos 


Po = ZP: + 3P, + 2P,. 


Geométricamente, esto significa que el vector Py “atraviesa” el 
paralelepipedo construido con los vectores P,, P,, P (fig. 17, b). 
Presentemos las mismas ecuaciones en la base P,, P3, P,. Para 
ello haremos una transformación idéntica. En la columna de 
la variante x¿ hay coeficientes positivos. Encontramos el min 


2.302 | a 
EE y Escogemos como elemento pivote al coeficiente de xz 
en la primera ecuación. Resolvemos esta ecuación con relación 
a xa y colocamos la expresión obtenida en las otras dos. Entonces 
las ecuaciones toman la forma: 


De las últimas ecuaciones se deduce que en la base P,, Py, P,, 
el vector Py se representa en la forma P, =0-P,+P,+P,. 
En este caso tiene lugar la degeneración ya que el vector Py se 
encuentra en el plano básico determinado por los vectores P, y P,. 
Por eso P, está en el lado del paralelepipedo formado por los 
vectores P,, P3, P, (fig. 17, c). Notemos que sí hubiésemos infringido 
las reglas de las transformaciones idénticas, hubiésemos llegado 
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a una base en la que P, no podria ser representado como una 

combinación lineal no negativa de vectores de base. Por ejemplo, 

si se toma por elemento pivote no al que corresponde el 
E ES j e 

min E 3 2) sino al coeficiente de la variable x, en al tercera 

línea, entonces después de las transformaciones, las ecuaciones 

tendrán la forma: 


Xa4 —_ 2 = (x3 + 2x5), 
x, = —2 — (—2x3 — 3x5), 
X23 = —3 — (-3x3 -— 7x5). 


De estas ecuaciones se deduce que en la base P,, P,, P,, el vector 
P, se representa como una combinación negativa que contiene 
coeficientes negativos 


P, = —2P, — 3P, + 2P,, 
y la solución resulta no positiva: 
Xx = —2; X2 = —-3; x3=0; => 2: xs =0. 


En este caso el vector P¿ no atraviesa el paralelepíipedo formado 
por los vectores P,, P,, P, (fig. 17, d). 
Así que la infracción de las reglas de las transformaciones 
idénticas definidas anteriormente puede llevar a una base negativa. 
- Veamos un ejemplo de un sistema de ecuaciones, las transfor- 
maciones idénticas del cual pueden llevar a la periodicidad: 


x1 =0— (X4 — X5 — X6 + 3x3), 
X> 20 (2x0 13 30 40) 


X3 = 1 — (Xa4 +5 3X6 — 8x»). 


Examinemos las transformaciones de este sistema con relación 
a la tercera ecuación. En ésta, el coeficiente de x¿ es positivo 


; ; 0. 0. 1 
e igual a uno. Entre las relaciones (2: n +) destacaremos 
la primera que es una de las menores entre ellas. Por eso se puede 
tomar como elemento pivote el coeficiente de la variable x¿ de 
la primera ecuación. Después de la resolución de la primera 


ecuación respecto a x4 y de la sustitución de x, por su expresión 
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en las otras dos ecuaciones obtendremos el sistema transformado: 
Xa4a =0— (xx, — xs + X6 + 3x3), 
x, =0- (2% + Xs + xo -— 5%) 
xx =1l-(-x, + 2x5 + 4x6 — 11x»). 


En la última ecuación en la parte derecha de nuevo hay coeficientes 
positivos de las variables, por ejemplo de xs. El elemento pivote 
en este caso es el coeficiente de x; de la segunda ecuación, puesto 


que el mín jo: + corresponde a la segunda ecuación. Efectuemos 


las transformaciones idénticas que corresponden al elemento pivote 
nuevo. Como resultado, obtendremos el sistema: 


Xs5 = 0 — (2% $ X2 + xs Pa 5%, ) 


10 (=x + X> +20) 
Xq = 1 — (—3x, — 2X» E Ne — X3). 


En la tercera ecuación aparece de nuevo un coeficiente positivo 
en la incógnita xs y se puede tomar como elemento pivote el 
coeficiente de la variable x¿ en la segunda ecuación. Después 
de la transformación idéntica tendremos: 

X6 —= 0 == (—2x,; + 2X» + 2Xa a 4x7), 

x5s=0-(. x1-2x2-3x4+ x7), 

X3 = 1 —( Sx1 ea Ax» e 2Xa4 + 3x3). 


Continuaremos la sucesión de estas transformaciones subrayando 
cada vez en la tercera ecuación uno de los coeficientes positivos 
y destacando con caracteres gruesos su respectivo elemento pivote. 
Obtendremos así el sistema inicial de ecuaciones expresado en 
diferentes bases propias positivas: 


x =0-—( xi — 2x, — 3X4 + Xs), 
Xp =0 — (2x, — 6x, — 10x4 + 4x5), 
Xx3=1-—(Qx, +2x, + 7x4 — 3x5). 


SO 


x¡ =0-— —3x, — Sxa + 2x5 ++) 


x =0-( X + 2x4 — xs 0) 


x3=1-—( 8x2+17x4-—7x5=  Xg). 


Xx3=1l-—(x4+xs+ 3x6 — 8x7). 


Según esto, después de seis transformaciones idénticas volvemos 
al sistema inicial de ecuaciones. 

Asi, la sucesión de las transformaciones idénticas en las bases 
P,P,P,>P,P,P3>P5P,¿P,—>P¿P;sP3>P,P¿P3—P,P,P3—P;,P),P, 
lleva a la periodicidad. 

No obstante, la existencia de la degeneración no conduce 
obligatoriamente a la periodicidad, lo que se deduce del ejemplo 
siguiente: 


Xx =0-—( x4—2x5+4x56+ x7—3xg + Xo), 
x =0-—(2x4-— xs5+ Xx — 2x7 — 2xg + 4x0), 
X3 =0—(—xXx4 + 2x5 — 3x6 + 3x7 —4xX3 + Xo), 
O=1-(3x4— Xx5—2x6+ X7— Xg+ Xo). 


En la última ecuación, el coeficiente de la variable x, es positivo, 
y como elemento pivote se puede tomar el coeficiente de x4 
en la primera ecuación. 

Efectuando las transformaciones idénticas tendremos: 


Xa =0—(. x1—2x54+ 4x3+ x7-—3xg+ Xo), 
x, =0— (-2x1,+3x5— 7x6 — 4x7 + 4xg + 2x0), 
x3=0-=( xi + Xó6 + 4x7 — 7Xg + 2x0), 
0 =1-— (—3x, + 5x5 — 14x5 — 2x7 + 8x3 — 2X0). 
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Continuando este proceso obtendremos: 


2 1 7 4 4 
Xg = == (a. + 3 X2 3% 7% + 3 X8 + 3 10) 
1 Z 2 5 1 7 
m0 (nta pt 30) 
x3 =0-—( Ni + x6+ 4x-— 7Tx3+ 2x9), 
1 5 7 14 4 16 
0 =1- — X1y =— X> —= —X6 FT X> F—Xg =——Xo l. 
E E E E A ) 
Y después 
1 
| q + a qe A 20) 
1 E 1 6 3 
xs =0-( E O 3 X3 3 X6 2 ve 3 1) 
1 Z 5 3 39 38 
=0-| —x1 +2 — X3 — —Xg H-— Xo ), 
de Te IR CNO VAGO EA ) 
A A: A O e o dp O 2, 
O E CALA y E AS 


En la última ecuación el coeficiente positivo de xg es simultánea- 
mente elemento pivote, de modo que no se ha formado una 
sucesión infinita de transformaciones idénticas. Esto muestra que 
la degeneración no ha llevado a la periodicidad. 

Efectuando las tramsformaciones idénticas que corresponden 


al elemento pivote axg == y admitiendo que el valor de las 


1 
variables no básicas es igual a cero, encontraremos una de las 
soluciones básicas del sistema de ecuaciones inicial: 

x1=0; x2=0; x3=0; x5¿=0; x9=0; 


13. “Xx AE Xx sais XxX = >— 
18” 5 38 > 7 , 8 . 


X4 


Es de señalar que la composición de ejemplos con periodicidad 
presenta grandes dificultades puesto que para la realización de 
la periodicidad, los coeficientes del sistema de ecuaciones lineales 
deben satisfacer un gran número de condiciones y solamente una 
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minuciosa selección de los coeficientes puede llevar a la formación 
de un ciclo. 
Mostremos que si se tienen sólo dos ecuaciones 


ñn 
x1 =0-— (a1axs + U14X4 + E 00%) 
¡j=5 


je 
(19) 


n 
X23 = b-— (dass, + l24X4 + ya 0290) 
153 


no puede haber periodicidad. La matriz de los vectores que corresponde 
al sistema (19) tiene la forma 


P, Pz Po Py P, P; 
1. 0 0. 413 44 4, j 
0 1 ba 424 4) 

Si se admite que b > 0, entonces los vectores P, y P, forman una base. 


El valor de los componentes del vector Pz en esta base se puede 
representar con las determinantes 


di3 = 


a13 O . am =|! 013 
d23 1 0 az 


Sean 4,3 Y 4,3 mayores que cero. Eso quiere decir que se puede someter 
la ecuación (19) a transformaciones, pasando de la base P,, P, a la base 
P,, P> y tomando el elemento a,3 como elemento pivote. En la nueva 
base Py, P,, los coeficientes de la ecuación (19) toman el valor determinado 
por las fórmulas (7) del capítulo anterior. Así, los valores de los 
coeficientes de la variable x¿ tendrán la forma: 


' E 1 
en la primera línea — | %14 li 
424 1 
, a a a 0 
en la segunda línea — | “9 “1% | donde A=| “*? !/ 
423 Az4 da3 1 


Si los numeradores de estas expresiones son positivos, entonces en la base 
P,, P,, los coeficientes de x¿ serán positivos. Puesto que la relación 
mínima en el sistema (19) de nuevo corresponde a la primera ecuación, 
tomando el coeficiente de x¿ de la primera ecuación como elemento 
pivote se puede repetir la transformación del sistema de ecuaciones 
pasando de la base P,, P, a la base P,, P,. 

Supongamos que después de cada transformación se puede realizar 
la siguiente, de tal manera que la sucesión de las bases 


P;: P, > P,, P,>P,, P),>..>P,, P,>P,, P, 
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forma un ciclo. Entonces, cada una de las determinantes, en el orden 
siguiente: 


(P,, P2) (P3, Pz) (Pa, P,) (Ph-1P2)  (P¡P2) (PP) | 
s 1 y so? s y (20) 
(PP3)  (P3Pa)  ... (P.-2Pi-1) (P,-1Pi) (P,P1) 
debe ser positivo (las determinantes de la fila inferior corresponden 
a la sucesión de los coeficientes positivos en la segunda ecuación, las 
determinantes de la fila superior, a la sucesión de elementos pivotes). 
En el caso dado los vectores P; son bidimensionales y por eso el valor 
de cada una de las determinantes en la sucesión (20) coincide en signo 
y magnitud con el respectivo producto vectorial Entonces la exigencia 
de que todas las determinantes (20) sean positivas es equivalente a la 
exigencia de que el vector Pz caiga entre los vectores P, y P, que 
forman un ángulo agudo; el vector P¿, entre los vectores Pz y P,, etc. 


Fig. 18 


Continuando estos razonamientos llegaremos a que el vector P, debe 
encontrarse entre los vectores P, y P, lo que no puede ser (fig. 18). 
La contradicción también puede ser fácilmente revelada en forma analítica. 
Escribamos el sistema de igualdades vectoriales: 


P, = a3¡P, + a32P,, 
P, = %43P3 + 9%42P, 
P; = ds¿P, + a52P,, (21) 


P. = %-1Px-1 + %%2P, 
P, = Pr + a,P,. 


Multiplicando vectorialmente cada una de las igualdades (21) una vez 
a la derecha por P, y la otra vez a la izquierda por el segundo vector 
que entra en la parte derecha de cada ecuación, nos convencemos de que 
todos los coeficientes de las igualdades (21) son positivos, puesto que, 
según la suposición, todos los productos vectoriales respectivos son también 
positivos. Así pues, para la ecuación Pz =au3¡P, + %3,P, tenemos que 
(P¿P,) = 93, (P,P>) + 9%, (P¿P,), por consiguiente %3, > 0. Multiplicando 
esta misma ecuación a la izquierda por P, obtendremos que «32 > 0. 
Después, en la expresión de P, sustituimos P,., por su expresión, en la 
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fórmula obtenida sustituimos P,-_, por su expresión, etc. Como resultado 
tendremos P, = Y1;P, + Y2;P>> Y1x > 0; Y2x > 0, de donde 
1 
A A 
Y1k Yik 


Comparando la igualdad obtenida con la última igualdad del 
ist 21 
da P, =01P, +012P>, 0%>0; 012 >0, 
llegamos a una contradicción. 

Mostraremos que para formar un ciclo, en el caso de cualquier 
número n de ecuaciones, n> 2, hay que efectuar no menos de seis 
transformaciones idénticas. Puesto que examinamos sólo el caso de una 
sustitución cada vez en las bases, se puede esperar el regreso a la base 
inicial después de 2k transformaciones, donde k=1, 2, 3, ... El regreso 
en el caso de k=1 es imposible, puesto que las determinantes 
correspondientes a este caso tienen la forma: 


(PyPz .. Pi... Pp); (PiPz .. Pj ... P,);5  (P¿Pz ... P; ... P,) 


s T y T 
(P,P, ... P; ... P); (P,P, ... P; ... P,). 


Las dos últimas determinantes se diferencian en una permutación 
de los vectores P, y P; por lo que tienen diferentes signos. Eso contradice 
la suposición de su positividad. 

Para demostrar la imposibilidad de la formación de un ciclo cuando 
k = 2 hay que mostrar que ninguno de los caminos de las transformaciones 
idénticas que se presentan más adelante es irrealizable. 


1. (P,P, .. P, .. P;, -. P,) >(P¡P .. P; .. P;, .. P,)> 
>(P¡P, .. P, .. P,, .. P,)=> 
>(P¡P, .. P, .. P;, .. P,) => 
>(P,Pz ... P; ..« P;, «. P,); 

2. (P¡P, ... P, .. P; .. P,) >(P¡P> .. P; .. P;, ... P,) > 
>(P,P .. P; .. Pj .. P,) > 
>(P,P, .. P, .. P;, .. P,) => 
LAS A A 

3. (P,P, 0. Pj .. P; . Pp) >(P¡Pz .. P; .. P;, .. P,)> 
>(P,P, .. P, .. P; .. P,)> 
>(P/Pz ... P; e Py. P,) > 
>(P,P, ... P; .. P,, .. P,). 


El primer camino de transición de base a base es imposible, puesto 
que se reduce al caso de las dos ecuaciones. 

El segundo camino tampoco es posible puesto que las determinantes 
primefa y última se diferencian en una permutación y por eso no pueden 
ser las dos positivas. 

Queda por demostrar la imposibilidad del tercer camino. Para eso 
hay que convencerse de que la suposición sobre la positividad de todas 
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las determinantes del r-ésimo orden anotadas a continuación lleva a una 
contradicción: 


l Sta Ly | sae 1, ta P 1 a 11 Pb 1 e.» PF 


| ] : a ] 
ti ó , A IA 
ly ij ' 8 E SETE 
; ; 
12 eE > : 4 in as 
F a 1 a a l 


Y, ] y é UB ? ri 


[> E 


la) hi ) 
Ari, : l aj | 
rj1 
1 ly 15) / 
| e, A, 
A, 12 ; 
Si 1 Ada 
Aj, Ari 
l Ey bh) r 
| A : ij 
Ml Gio y 5 
: , 
ba a <= ue 
| | 
á Aj, Ar; 
| 14 l , 
[ % 
| M1, 
! NA 
Aj, Í 
a 


S6 


Resolviendo estas siete determinantes tendremos: 


1) A > O; 4) Ar ¡ Api, — Aj jj, > O; 
2) 4y;, > 0; 5) aj, > 0; 
3) Mi, Mi,j, — Gi, %,y, > 0; 6) 9, j,Ayj, — Ari), < 0; 


7) 4), <0. 


De 1), 2), 4) y 7) se deduce que a, , < 0. De 2), 5), 6) y 7) se deduce que 
a. <0 


Introduciremos las designaciones: 


zi 


A, es Ar, A 
Ari, = ke, Aia NA k,n» 
Aja PS Ap AA qí18. 


€, K, q, €.» Kk1, q1, son números positivos. 
Entonces de la 3) tenemos «qe — k,qyn? > 0 | qe? — k,q,n? > 0, 
de la 4) tenemos — en + k,ken > 0 | kk, — 1>0, 
de la 6) tenemos qn? — kqe? >0 | qn? — kqe? > 0, 
qe” > k,gin”, 
qun? > kqe” 
qq1e.m? > k,kq,qe?n?,  1>kik. 
Comparando la desigualdad 1>k,k con la obtenida anteriormente 
k,k > 1 observamos una contradicción. 


Así que la periodicidad puede ocurrir solamente después de que 
se efectúen no menos de seis transformaciones. 


$ 8. METODO DE CÓMPUTO DE UNA SOLUCIÓN 
NO NEGATIVA DE UN SISTEMA DE ECUACIONES 
ALGEBRAICAS LINEALES 


Sea dado un sistema de m ecuaciones algebraicas lineales con 
n variables: 


411X1 +019X2 +... + a1 5%; A A bi, 
A71X1 + 422X2 + ... =P 4 +... + 029Xpn = b,, 


AmiX1 + Am2X2 + 0. + AmpXj + 2 + Amon Xp = Dip > 
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Sin limitar la generalidad podemos considerar que b,> 0 (esto 
se puede lograr siempre multiplicando las dos partes de las 
ecuaciones respectivas por —1). 

Presentemos el sistema (13) en la forma 


0 — b; == ($ 09%) Í — 1, Ze ss. IM. (22) 


j=1 

Si en el sistema (22) hay alguna variable que entre sólo en una 
ecuación y el coeficiente de la cual, dentro de los paréntesis, tiene 
el signo «+», entonces se podrá resolver esta ecuación con relación 
a esa variable. 

Supongamos que las ecuaciones del sistema (22) se han resuelto 
con relación a todas las variables tales. Entonces, después de la 
respectiva numeración de las variables, el sistema puede ser escrito 
en la forma: 


Xx1 = bj — (A17, 41%1, +1 * O9, 42%, 42 4 o» + QyXj +... + yg Xp), 
X2 = b, — (Aa, +1%1,+1 + Ga1,42%1,42 4 0 E QyX¡ + o... + ay Xp), 


.. . . ...o.. . ..oo.o.ox2.02.02..bnxÁxoxoxÁr£Ár£202e..o.xP202-2.2..2x2b0..o. o . . . .—«—2-...0.. 0... ....- 0.0..202Pb04.-.—..o 


Xt, = Pro — (Moto + 1941 4 yy) 2X19 42 7 0 Oy Em E Xp), 
A AA AO A 
... + Ar, +1, kXk), 
0= Bm = (Am, 1,91%) +1 4 Gm 1, +21, 42 E 00 E OmppXj E a. E AX) 


(23) 


donde 
pa 1, Ze ..., lo; y = 1, 2; ..., YO> 


l+Yo=m; lo+k=n; b,>0; b,>0 


A cualquier ecuación del sistema (23) que no esté resuelta con 
relación a alguna de las variables la llamaremos 0-ecuación. 
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Asi que cualquier sistema de ecuaciones algebraicas lineales 
puede ser reducido a la forma (23)*”. 

Para encontrar una solución no negativa del sistema (23) lo 
someteremos a unas transformaciones idénticas que satisfagan las 
condiciones siguientes: 

1. Hallamos la O-ecuación, en la que el miembro libre b, es 
mayor que cero. Si no hay una 0-ecuación tal, entonces los valores 
de las variables 


Xx] = by; X= 0 (1 = E 2, Po lo; j= lo + 1, ds k) 


forman una solución no negativa del sistema (23). Supongamos 
que ésta sea la ¡-ésima ecuación. 

2. Señalamos en la ¡-ésima ecuación el coeficiente positivo a; e.) 

3. Hallamos el elemento pivote a, y efectuamos en el sistema 
(23) las transformaciones idénticas correspondientes. 

4. Empleamos la ¡-ésima O-ecuación para las transformaciones 
siguientes del sistema hasta que lo resolvamos o establezcamos que 
el sistema (23) es incompatible **. 

5. Después de resolver la ¡-ésima 0O-ecuación hallamos la 
0O-ecuación siguiente en la que el miembro libre sea positivo 
y efectuamos con ella operaciones análogas. 

6. Continuamos este proceso hasta que nos libremos de todas 
las O-ecuaciones. 


Observaciones: 1. El número de ecuaciones en el sistema transformado 
puede resultar menor que m, puesto que al introducir en las demás 
ecuaciones la expresión de la variable resuelta algunas de estas ecuaciones 
se pueden convertir en identidades O =0, las que habrá que excluir del 
tratamiento ulterior. 

II. Supongamos que en el sistema inicial (23) o en el sistema 
transformado se ha encontrado una ecuación resuelta con relación 
a pa variable, la que tiene el miembro libre igual a cero: x;= 


— OB, Xo), y es tal que en la parte derecha todos sus coeficientes 
de je variables son no negativos. Como el papel de tales variables puede 
ser solamente de cero, para que las x, sean no negativas hay que 
igualarlas a cero en todas las ecuaciones con lo que el sistema examinado 


*) Si en el sistema (23) no hay ni una ecuación 
resuelta con relación a alguna variable, entonces ly =0 y el sistema (23) 
está compuesto sólo de O-ecuaciones. 

**) Si en la ¡-ésima ecuación, las incógnitas no tienen 
coeficientes positivos, el sistema (23) resulta as: ya que la 
ecuación 0 = b; — =(2 aX), en la que b,> 0 y todos los a;; < 0, no puede 
ser satisfecha con ningún valor no negativo de las variables x;,> 0. 
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se simplificará —considerablemente. En forma similar excluimos las 
O-ecuaciones que tienen el miembro libre igual a cero y los coeficientes 
diferentes de cero de las variables en la parte derecha de igual signo. 
Tomamos el valor de estas variables en las demás ecuaciones igual a cero. 


Veamos a qué resultados pueden llevar las transformaciones 
idénticas consecutivas que satisfacen las condiciones 1-6. 

a) Puede ocurrir que después de un número finito de trans- 
formaciones idénticas el sistema se libre de las 0-ecuaciones. 
Evidentemente, el sistema (23) en este caso es compatible. 
El conjunto de valores de las variables obtenidas al igualar a las 
variables no básicas a cero y a las básicas a los miembros libres 
en el sistema que no contiene O-ecuaciones, es una solución 
no negativa. 

b) Puede ocurrir que después de un número finito de trans- 
formaciones idénticas se ponga de manifiesto que la 0O-ecuación 
empleada se convierte en una ecuación de forma 


0 = b; — 03 Aj¡X 5), 


donde b;>0, a;;¡<0 para todas las j. En este caso el sistema 
es incompatible. 

c) Por último, puede ocurrir que el sistema no se libre por 
completo de las O-ecuaciones, y las condiciones del empleo de las 
transformaciones idénticas no se alteren. Como al emplear 
consecutivamente las transformaciones idénticas, el número de 
O-ecuaciones no puede aumentar, en este caso, cierta 0O-ecuación 
posee la propiedad de que ella siempre tendrá por lo menos 
un coeficiente positivo de las variables en la parte derecha pero 
nunca le pertenecerá el elemento pivote. Cambiando esta O-ecuación 

0 =b;, — (Dax) 
por la ecuación 
b; 1 
y o > E AjpjX5), 
donde e es un número positivo tan pequeño como se quiera, y, 
estudiando la última ecuación junto con las resueltas veremos que 
en el caso c) ocurre la periodicidad 

Lo dicho anteriormente da la posibilidad de hacer la siguiente 
afirmación: 

Sea cual fuere el sistema de m ecuaciones algebraicas con n 
variables, el empleo consecutivo de transformaciones idénticas permite 
después de un número finito de pasos (excluyendo los casos de 
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periodicidad) determinar si el sistema en la región de los valores 
no negativos de las variables es compatible. En el caso de la 
compatibilidad del sistema, en el último paso de las transformaciones 
se forma cierta solución no negativa. 

Una importante propiedad del método presentado es la 
posibilidad de su empleo para cualquier sistema de ecuaciones 
algebraicas lineales, incluyendo también el sistema linealmente 
dependiente, sin que haga falta separar de antemano las ecuaciones 
linealmente independientes. 

Además, como cualquier sistema de desigualdades lineales al 
añadirle variables no negativas puede ser reducido a un sistema 
de igualdades, este método también es utilizable tanto para 
determinar la compatibilidad de sistemas de desigualdades, como 
para encontrar una de las soluciones no negativas en el caso 
de su compatibilidad. Aquí tampoco es necesario separar al poliedro 
de las soluciones en su “forma pura”. 

En la fig. 19 se representa el esquema de bloques del programa 
para calcular la solución no negativa de un sistema cualquiera 
de ecuaciones algebraicas con ayuda de computadoras electrónicas. 

En la computadora se introducen los coeficientes de las 
variables en el sistema inicial y los miembros libres. Además, 
al resolver problemas de programación lineal, se introducen los 
coeficientes de la forma lineal que se minimiza. La computadora, 
a base del programa, efectúa la sucesión de transformaciones 
idénticas y comprueba las condiciones de compatibilidad del sistema. 

En el caso de incompatibilidad del sistema, la computadora 
da una señal indicadora y detiene su funcionamiento. Si el sistema 
es compatible, la computadora calcula la solución no negativa, 
la imprime (si es necesario) y pasa al cálculo de la solución 
Óptima. 

Veamos algunos ejemplos. 

Ejemplo. Determinar si es compatible el sistema en la región 
de los valores no negativos de las variables: 


1-x1 + 3x2 + 3x3 =2, 
l-x, + 2x7 + 4x3 =7, 
—3x,+3x2+ x3=3. 


0=2-—(1-x, + 3x, + 3x3), 
0=7-—(1-x, + 2x, + 4x3), 
0=3 — (-3x, + 3x, + X3). 


el sistema es 
incompatible 


Stop (parada del 
funcionamiento 
de la máquina) 


Paso al cálculo de la 
solución óptima 


Salida de la solución 
no negativa 


Se comprueba la existencia de coeficientes 
positivos de las incógnitas en las ecuacio- 


nes del sistema 


Elección del primer coeficiente positivo 
4, ¡+ €n la primera 0-ecuación 


Elección de los coeficientes positivos a; 
en la columna j* entre todas las ecuacio- 
nes del sistema 


, . b; 
Cálculo de las relaciones ——— 
Aij* 


bis 


Búsqueda del valor mín entre los 


Ajxj* 


números 
Ap 


Resolución de la ecuación ¿* con relación 
a la variable Xp 


Traspaso de la ecuación resuelta al primer 
lugar del sistema de ecuaciones 


Registro en la memoría del número ¡* 
de la variable xj: 


Exclusión de la variable Xy» de todas las 
ecuaciones del sistema 


Verificación de la existencia en el sistema 
de identidades tipo 0 =0 y reducción 
correspondiente del número de ecuaciones 


Verificación de la existencia de 0O-ecua- 
ciones 


Fig. 19 
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En la primera 0-ecuación, el coeficiente positivo de x,, igual 
a uno, es simultáneamente el elemento pivote: 


x, =2 — (3x, + 3x3), 
O=5 — (—x> + X3), 
0=9— (12x, + 8x3). 


En la segunda 0O-ecuación el coeficiente de xz es positivo; 
el elemento pivote pertenece a la primera ecuación: 


El sistema es incompatible, puesto que en la última ecuación 
el miembro libre es positivo y los dos coeficientes dentro del 
paréntesis, negativos; por eso, a la segunda ecuación no le satisface 
ni un punto para el que x, > 0; x,>0. 

Ejemplo. Hallar una solución no negativa del sistema de 
ecuaciones: 


0 =3 — (2x1 + x, + 2x3 + 2x4), 
(1) 0 =3 — (-2x, + x, + 4x3 + 2x4), 
0 = 6 35 (—2x; +2x» + 7x3 + Axa). 


Efectuamos las transformaciones idénticas del sistema: 


X1 2 (a + X3 +) 


(2) 0=6-— (2x, + 6x3 + 4x4), 
0 =9 — (3x, + 9x3 + 6x4). 


Una de las soluciones no negativas del sistema es 


Xi =—3 *Xxai=1l; x,=0; x¿4=0. 
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El método de cálculo de la solución no negativa de un sistema 
de ecuaciones explicado anteriormente puede servir como uno de los 
procedimientos de definición del rango de una matriz. Efectivamente, sea 
dada una matriz 


411 12 +++ Aj, 
41 422 +++ 42, 
Amt Am2 ++ Am 
Consideremos sus columnas como los vectores P,, P,, ..., P, y escribamos 


la ecuación 


0 = Po (2 «p,) 


j=1 


en la que haremos P, = P,. Es evidente que esta ecuación tiene solución 
no negativa (por ejemplo, x, =1; x¡=0; j¡=2, 3, .., n). Representando 


n 
la ecuación vectorial 0 = P, — ( ye «P,) en forma de un sistema de m 
j=1 


ecuaciones y sometiéndola a transformaciones idénticas nos libraremos 
de las O-ecuaciones. El número de ecuaciones resueltas después de la 
exclusión de las O-ecuaciones determina el rango de la matriz. Como se 
puede demostrar, el número de transformaciones que es necesario efectuar 
también es igual al rango de la matriz. 

Ejemplo. Hallar el rango de la matriz 


2 —4 3 1.0 
lt —2 l -4 2 
0 LL 4 3.0 
4 7 4 -4 5 


Efectuaremos la sucesión de transformaciones identicas: 


0=2-—(2x, —4x,+3x3+ Xa), 

O=1-—(x,¡-—2x2+ X3 + 4x4 + 2x5), 
0=0-( X) — X3+3X4+ Xs), 
O=4-— (4x, — 7x2 + 4x3 — 4x4 + 5x5) 


xy =1—(—2x2+x3-— 4x4 + 2x5), 
0=0-( x3 + 9x4 — 4x5) 
0=0-( X2=X3+ 3x4+ Xs). 
0=0-( X + 12x4 — 3x5) 


1*" paso 
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x3=0-( 9x4 — 4x5). 

x¡ =1-—(—2x,-— 13x4 + 6x5), 2% paso 
0O=0-(  x>+ 12x4 — 3x9), 
0=0-(.  x2+ 12x4 — 3x5) 

x, =0-— (12x4 — 3x5), 

x3=0-—( 9x4 — 4x9), 3 paso 


Xx 1 = l—1ixz, 
0=0 


Rango de la matriz r = 3. 


$ 9. RESOLUCIÓN DEL PROBLEMA 
DE PROGRAMACIÓN LINEAL 


El problema de programación lineal, como fue dicho anterior- 
mente, puede ser planteado de la siguiente manera: 
se da un sistema de m ecuaciones lineales con n variables 


411X1 + 012X2 +. + Oy¡Xj 4... + OinXp =D, 
431X1 + 022X2 + «+ 0a¡X¡ + ... + 025Xp =D, 
Aj1X1 + QAj2X2 + .. as aX; +... + AinXn = b;, 


AmiX1 + Am2X2 +... + AamjiX; +... + UmnXn = Dr 


y la forma lineal de esas mismas variables f =C,x, + C3X) + ... 
+ CjXj + ... + CpXn. Hay que hallar entre todas las posibles soluciones 
no negativas 


(Mir a) 


del sistema (13) una solución (xi; Xx3; ...3 x,) con la cual f 
adquiera el menor valor posible. 

Para resolver este problema calcularemos una de las soluciones 
no negativas del sistema (13) con el procedimiento explicado 
en el parágrafo anterior. Como resultado, el sistema (13) se 
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transformará en un sistema equivalente *): 
? , , 1 ,? !? , 

Xx =0D1 — (A1r+1Xp+1 + Mir 2Xp +2 4 0 E AX... + Xp), 
ro pr / ! , / o, ! 

x2 = Da — (0774 1Xp41 + Mary 2Xp42 4 0 E ¿GE o. + aX), 


(24) 


_,/ / / / Ú Ps Im, 
Xy = b, NN (Ap, + 1Xp+41 + Arr+ 2 Xp 42 700 + a, X; ho + Ark Xp), 
i=1, 2... ri r<m; r+ks<mn. 


Sustituyendo la expresión de las variables básicas en la forma 
lineal e introduciendo las designaciones respectivas para las 
magnitudes constantes, presentaremos f en forma de una función 
lineal 


k 
a — Ñ A 
f=F (E, 6%) (25) 

En concordancia con el método simplex, para hallar la solución 
óptima sometemos el sistema de ecuaciones (24) y la función 
lineal (25) a transformaciones idénticas hasta: que se infrinjan las 
condiciones de su realización, o sea, hasta que desaparezcan los 
coeficientes positivos de las variables en la forma (25). Pero como 
lo expuesto en los dos primeros parágrafos de este capitulo sobre 
la periodicidad se refiere también al problema de programación 
lineal, se puede afirmar lo siguiente: 

El empleo de un número finito de transformaciones idénticas 
permite (excluyendo los casos de periodicidad ) encontrar una solución 
no negativa del sistema de ecuaciones con la cual la forma toma 
el valor minimo. 

Se han propuesto varios métodos con el fin de evitar la 
periodicidad. Algunos de estos métodos al efectuar las transforma- 
ciones idénticas, exigen la realización de operaciones complemen- 
tarias cada vez que ocurre la degeneración. Cuando ocurre 
la degeneración, el elemento pivote no se determina en forma 
única. Las operaciones complementarias, en este caso, en cada 
paso de las transformaciones idénticas, permiten elegir el elemento 
pivote de tal manera que la periodicidad se excluya. El caso 
de la degeneración, como ya se ha indicado, siempre tiene lugar 
cuando el vector P, pertenece a uno o a varios hiperplanos básicos. 


*) Después de una numeración correspondiente de las 
incognitas. 
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El análisis de la resolución práctica de problemas de 
programación lineal confirma la suposición sobre la pequeña 
probabilidad de la periodicidad, además la composición de ejemplos 
con periodicidad ha encontrado dificultades bastante considerables 
debido a la necesidad de satisfacer un gran número de condiciones 
complementarias. 

Aquellos pocos casos en los que se forma un ciclo pueden ser 
hallados directamente al resolver el problema en forma manual 
(sin computadora). 

Por ejemplo, planteando el problema del cálculo de una solución 
no negativa del sistema 


x1 =0 — (x3 — X4 — Xs + 3x6), 


X2 20 (2x9 +1), 


en el que la función lineal 
f =3—(x3 + xa + 3x5 — 8xg) 


alcanza su menor valor, y efectuando una sucesión de trans- 
formaciones idénticas como en el ejemplo examinado en pág. 49, 
observamos que hemos llegado al grupo inicial de las variables 
básicas. S1 no se presta atención a esto, el proceso de resolución 
puede resultar interminable, al mismo tiempo el valor de la 
función lineal quedará siempre constante. 

Al observar la periodicidad hay que cambiar el orden de las 
transformaciones idénticas eligiendo otro elemento pivote. Así, en 
el ejemplo examinado, procediendo tal como se explica a continua- 
ción, establecemos que la función lineal f no tiene límite inferior: 


X1 = 0 — (x3 — X4 — X5 + 3x6), 
x =0-— (2 — Xa xs + 10) 


f=3—(x3 + x4 + 3x5 — 8x6), 


1 1 1 
x =0- (q y *4— reia > 12) 
1 1 3 5 
y 0 (09 q +0) 
1 3 13 17 
=3=| -—x Xa FF Xs == 
/ > a > 4 + E %,) 
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En el caso dado, el valor de la forma puede ser disminuido 
a cuenta del incremento del valor de x,, pero como en la columna 
de los coeficientes de la variable x, no hay coeficientes positivos, 
O sea, que no hay ni un elemento pivote, esto quiere decir que 
no existe ningún límite para aumentar la variable x,. Así la 
función f se hace ilimitada en su parte inferior. 

Al resolver estos problemas en las computadoras electrónicas 
hay que prever en el programa la posibilidad de descubrir el 
regreso a grupos de variables básicas que ya se habían encontrado. 
Sólo en el caso de que se revele la periodicidad, se debe utilizar 
la parte del programa en la que se prevén las operaciones 
complementarias que excluyen la posibilidad de esta periodicidad. 

La operación complementaria más sencilla es la elección de otro 
elemento pivote en el paso en el que se ha observado la repetición. 

Ejemplo. Hallar el menor valor de la función lineal 


f =5x, — 10x, + 7x3 — 3x4 
en el conjunto de soluciones no negativas del sistema de ecuaciones: 


0 => (15 + x2 + 7x5 + 2x4), 


0 => — (2x1 — xa + 3x3 + 3xa), 


0O=4-— (2x1 + 2x, + 8x3 + Xx4) 


X2 =2= (1 + 4x3 +42) 
3 
— (2%, + 3x3 + 21) 


BES 

2 

7 7 
0==-|-— 7x3 + X 

> X1 + oe > .) 


a En 
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1 1 
X3 ES 6 x), x4 =1 — (2x3), 


3 3 
Xx =0 — (3%), X2 A — (3x3), 
3 


f=> (Pa) f = 18 — (-43x3) 


Es —18 con X1 =0( 


x3=0 
Xx =1 
PEN 
2 


Ejemplo. Entre las soluciones no negativas del sistema de 
desigualdades: 


—20x, + 12x, — 15x3 < 60, 
x1 + 2x, — 3x3 < 6, 
3x1 + 6x2 + 4x3 < 12, 
—20x, — 15x, + 3x3 < 60, 
10x, + 5x, — 2x3 > 10, 
6x, + 7x, + 42x3 > 42, 


hallar la solución con la cual la forma lineal f =x, +, +; 
alcanza su valor mínimo. 


Resolución. Introduciendo variables complementarias, conver- 
timos el sistema de desigualdades en el sistema de ecuaciones: 


y, = 60 — (-20x, + 12x, — 15x3), 
ya = 6 — (x1 + 2x7 — 3x3), 

ya = 12 — (3x, + 6x, + 4x3), 

Ya = 60 — (—20x, — 15x, + 3x3), 
O = 10 — (10x, + 5x2 — 2x3 — Js), 
O =42 — (6x, + 7x2 + 42x3 — y6) 
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y lo sometemos a transformaciones idénticas: 


7 6 42 
y 9 (Po +>3x 42 ya) 
yy =8-— La, +2x +10) 
ea dE: za -—=x, +21) 
0=12= —x, +12 Yo =3) 


Xa sii ESA ds 
3 TS 
pm 953 (1383, 665, _ 125 
s 6 54 2 "Tp?" 7?) 
2203 E DAS 
| 220 17 161 
a 17 161 
112 54 2 +379Ys* 752 ya), 
80 (265, 188) _ 35 
Ya 6 54 2779 257 7519 +6) 


Una de las soluciones no negativas del sistema de desigualdades 
es la solución x; = Xx =0; x3 =—. Expresemos ahora la 
forma lineal por medio de las variables secundarias: f = 
is E 


Como se puede advertir en la última expresión de la forma 
lineal, en la región de los valores no negativos de las variables 
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es ya imposible continuar disminuyéndola. Por eso la primera 
solución xy == x2 =0; x3 == ha resultado la óptima y el valor 


correspondiente de la forma lineal es f ,¡, = 2. 
$ 10. SOBRE UN PROBLEMA DE MIN-MAX 


Sean dados n vectores m-dimensionales P;, y el conjunto 
T = (t,j de los números reales £;, j = 1, 2, ..., n, que les corresponden. 

Examinemos el conjunto Y = (y) de todas las soluciones posibles 
no negativas y =(x;) de la ecuación 


x¡P, Ej xP, +. + xP, = Pos 


donde P, es el vector dado en el espacio m-dimensional que no 
es igual a cero. 

Convendremos en designar con x; a las variables x;x+ 0 en 
cualquiera de las soluciones y, y a los valores correspondientes 
b;, con b;. 

Con (t;,y designaremos el sistema de números t; que correspon- 
den a las variables x;+ 0 en la solución y. 

Acordaremos designar con t, al número mayor del sistema de 
números (f;,;. 

Así tendremos que 


L, = máx Litas > 
J 
Presentemos ahora el planteamiento del problema: 


Entre todas las soluciones posibles ye Y hay que hallar una 
solución y, para la cual el correspondiente 


opt > pe [máx tl js] > 
o bien 
opt => Ser (2,). 


Llamaremos óptima y designaremos con ygp; a la solución y 
(si ella existe) que satisfaga esta condición. Estos problemas se 
llaman problemas de mín-máx. (Un caso particular de este problema 
es el de transporte por el criterio del tiempo, que se estudia 
detalladamente en el capítulo IV). 
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S1 un SIsicma de m ecuaciones algebraicas lineales con n variables 


js bi, 1<i<m; 1<j<n, 
j=1 
n 
que corresponde a la ecuación vectorial Y x,P,= P,, se somete 
je 
a transformaciones idénticas, podremos después de un número 
finito de pasos obtener una solución óptima. 

Para ello es necesario: 

1. Calcular la solución inicial no negativa del sistema de 
ecuaciones. 

2. Hallar entre las variables básicas x, aquella que corresponde 
al mayor valor de t;, igual a £páx, 

3. Excluir del examen posterior las variables secundarias x;j, 
para las que t;¡> tpáx, 

4, Encontrar en la ¡-ésima línea, dentro del signo Y un 
coeficiente positivo a¡» y determinar cuál es el elemento pivote 
Aix, efectuando después las transformaciones idénticas. 

5. Repitiendo (si eso es necesario) la etapa 4 varias veces, 
sacar la variable x;, del conjunto de las variables básicas y excluirla 
del examen posterior. 

6. Repetir las etapas 2-5 hasta que, en cierto paso, en la 
línea correspondiente del sistema transformado de ecuaciones, todos 
los coeficientes de las variables se hagan no positivos. 

La solución básica obtenida en este paso es la óptima. 

En efecto, supongamos que después de cierto número de 
transformaciones idénticas, el sistema inicial toma la forma: 


E bs, —_ Q A; X5), 


.. . . .*.o.o...o.. .. . .o.o2.k 


Supongamos que la variable x, corresponde al mayor valor de 
tmáx. Entonces, si b, <0, excluyendo x, del examen posterior, 


obtendremos un sistema incompatible. (Si b, =0, de acuerdo con 


la observación II de la pág. 58 hay que tomar X;, =0 y también 


x¡=0 para las j para las que a;,,>0. Después se continúan 
las transformaciones idénticas). 
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CAPITULO In 


RESOLUCIÓN DEL PROBLEMA DE TRANSPORTE 
POR EL CRITERIO DEL COSTO 


A continuación se estudia uno de los problemas típicos de 
la programación lineal, el llamado problema de transporte. 
Al planificar la transportación de mercancias con frecuencia surgen 
cuestiones sobre la organización más racional de los transportes. 
En algunas ocasiones eso significa la realización del cálculo de 
un plan de transportes según el cual los gastos sean los minimos. 
En otros casos es más importante ganar tiempo y por eso se 
plantea el problema de la elección entre los posibles planes 
de transporte aquel que al realizarlo permita hacer llegar las 
mercancías a su usuario en menos tiempo. 

El primer problema ha recibido la denominación de problema 
de transporte por el criterio del costo: el segundo, problema 
de transporte por el criterio de tiempo. 

El primer problema es un caso particular del problema de 
programación lineal y puede ser resuelto con el método simplex 
el que se ha expuesto en el capítulo anterior. No obstante, las 
particularidades de este problema hacen que se resuelva más 
fácilmente con el método combinatorio. Cuando el número de 
puntos de partida y de destino es pequeño, éste método permite 
resolver semejantes problemas sin emplear computadoras. Cuando 
el número de puntos de partida y de destino es grande, el 
problema solamente puede ser resuelto con el empleo de las 
computadoras electrónicas. Por ejemplo, un problema de transporte 
con 30 puntos de partida y 40 puntos de destino se resuelve 
en la máquina “Strelá” en 25-30 minutos. 

En el capítulo III se expone el método combinatorio de 
resolución del problema de transporte por el criterio de costo 
y se presenta un esquema de bloques simplificado para resolverlo 
en las máquinas computadoras. 


$ 11. PLANTEAMIENTO DEL PROBLEMA 


Se puede formular el problema de transporte por el criterio 
del costo de la siguiente manera: 

Designemos con a;,, 4», ..., Am la cantidad de unidades de 
mercancías de la misma clase que se encuentran en cada uno 
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de los m puntos de partida, y con b,, b,, .., b, la cantidad 
de unidades de mercancías en cada uno de los n puntos de destino. 

Sea x;; la cantidad de unidades de mercancías que se planifica 
transportar desde el ¡-ésimo punto de partida hasta el ¡-ésimo 
punto de destino, y C;;, el costo del transporte de una unidad 
de mercancías desde el ¡-ésimo punto de partida hasta el ¡-ésimo 
punto de destino. 

Sea la cantidad de mercancias que se envía de todos los m 
puntos igual a la cantidad de mercancías que se necesita en los n 
puntos de destino. En tal caso se debe cumplir la condición 


uN (26) 
i=1 j=1 


Anotaremos las condiciones de problemas de este tipo en forma 
de tabla 2. Llamaremos solución del problema de transporte 


Tabla 2 


a la matriz X =(x;j) de un orden de mx mn, cuyos elementos 
no son negativos x;;> 0, que satisface las condiciones 


x=a (i=1, 2, ... m), (26) 
J 


3 xij=bj, (¡=1, 2, .., n). (26”) 


La condición (26') significa que del ¡-ésimo punto de partida se 
saca toda la mercancia; la condición (26”) significa que la necesidad 
del j-ésimo punto es satisfecha por completo. El problema se 


74 


reduce al cálculo de los valores no negativos de x;; 
el costo general de transporte 


m n 
Se 
i=1 


j=1 


¡¡» con los que 


sea el menor. 

Empleando la definición del producto escalar de las matrices 
se puede formular de otra manera el problema de transporte por 
el criterio del costo. 

Sea dada la matriz C= (c;;), donde c;, son números no 
negativos reales. Hay que hallar entre todas las soluciones X 
anotadas en forma matriz y que satisfagan las condiciones (26') 
y (26”) una solución con la que el producto escalar (CX) alcance 
el menor valor. Llamaremos óptima a la solución presentada en 
forma de matriz X que satisface esta condición. 


$ 12. SOLUCIONES BÁSICAS DEL PROBLEMA 
DE TRANSPORTE POR EL CRITERIO DEL COSTO 


Examinemos la matriz dada C =(c;)) y una matriz-solución 
arbitraria X = (x¡¡) (i = 1, 2, ..., m); (¡= 1, 2, ..., n). Consideraremos 
que m< n, puesto que en el caso contrario se las podría cambiar 
de lugar. 

Introduzcamos algunas definiciones. 

Llamaremos celda a un par de números reales ij, conjunto 
a un grupo cualquiera de celdas. A una sucesión de celdas en 
forma i;j,, ja, iaja, l2j3, .., la llamaremos cadena. Una cadena 
será cerrada si tiene la forma 

lijr, dida) daja> «e» Tio 1d 
Llamaremos ciclo a cualquier cadena cerrada. Un conjunto se 
llamará cíclico si contiene por lo menos un ciclo, y acíclico si 
no contiene ninguno. 

A cada celda ij le corresponde un y sólo un elemento x;; 
de la matriz-solución X, y también un y sólo un elemento c;; 
de la matriz de costos C. Por eso cualquier conjunto de celdas 
es simultáneamente un conjunto de los respectivos elementos, 
tanto de los x;; como de los c;;. 

Numeremos los elementos de la cadena cerrada O, yendo por 
ella, por ejemplo, en la dirección de las agujas del reloj. Diremos 
que los elementos de esta cadena que tiene números nones forman 
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la semicadena impar O", y los elementos con números pares, 
la semicadena par OP. 
Acordaremos designar la suma de los elementos c;; de, la 


semicadena impar con Y c;,;, y a la suma de los elementos de la 
o” 


ij> 


semicadena par, con Y C;; 
O? 

TEOREMA. Cualquiera que sea la solución X con un conjunto 
ciclico de elementos x;; diferentes de cero, existe una solución Y 
con un conjunto acíclico de elementos y;, tal que (CY) <(CX) 
y en la que el número de elementos de la solución diferentes de cero 
es menor que el número de tales elementos de la solución X. 

DEMOSTRACIÓN. Supongamos que los elementos x;; diferentes 
de cero de la matriz X forman una cadena cerrada 0O,. Comparamos 


las sumas Y) c;; y ) c;;; una de las sumas, por ejemplo YC; de- 
ez e; er 
berá ser tal que »c¡¡< ) c;;. 
er e; 


_ Formamos a partir de la matriz X =(x;¡) una nueva matriz 
= (x¡¡) cambiando elementos en la cadena O, de la siguiente 
forma: 


/ pa e ! a 
A E o o TA 


aja aga TA o iy = Xi, 
Aquí iij1, daja, ..., 1 j, son las celdas de la semicadena impar O”; 
iz, l2j3, .., 1f,, las celdas de la semicadena par O”; xpín es el 
elemento menor en la semicadena impar O”. Para todas las demás 
celdas ij tendremos x;;= X;;- 

Se puede advertir fácilmente que la matriz X”, obtenida a partir 
de la matriz X al transladar por la cadena O el elemento xpín, 
también es una solución del problema de transporte. El número 
de todos los elementos diferentes de cero de la matriz X” se hará 
menor por lo menos en uno. 

Como 


(CX) =(CX) + 0, Cij — 2. y, 
O? 


entonces, (CX” < (CX). 

Continuando la construcción de tales soluciones después de 
un número finito de pasos llegaremos a la solución Y, la que no 
tiene cadenas cerradas O de los elementos x;, diferentes de cero, 
y (CY) <(CX). El teorema queda demostrado. 
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COLORARIO. Para encontrar aunque sea sólo una solución óptima 
es suficiente investigar todas las posibles soluciones X con conjuntos 
aciclicos de elementos x;; diferentes de cero. 

Cualquier solución X, presentada en forma de matriz, cuyos 
elementos x;; diferentes de cero forman un conjunto acíclico, se 
llama básica. 

En general, en una matriz X que represente una solución 
cualquiera, puede haber distinto número de elementos x;, diferentes 
de cero. No obstante, independientemente de la solución X que 
tomemos, ella tendrá no menos de n elementos diferentes de cero, 
lo que se deduce de la condición (26”) del $ 11. Por otra parte, 
se cumplirá el teorema siguiente: 

TEOREMA. El número N de elementos  x;;, diferentes de cero, 
de cualquier solución básica satisface la condición n < N <m=+mn-— 1, 
en la que m es la cantidad de puntos de partida, y n, la cantidad 
de puntos de destino. 


DEMOSTRACIÓN. Previamente nos convenceremos de que son 
justos dos lemas. Construiremos en un espacio (m + n)-dimensional 
m Xx n vectores y los pondremos en una concordancia mutua de un 
significado único con el conjunto de todas las celdas ij de la 
matriz C =(c;;j). Pondremos la celda ij en concordancia con el 
vector P;; que tiene todos los componentes iguales a cero excluyendo 
el ¡-ésimo y el m-+j-ésimo, cada uno de los cuales es igual 
a la unidad. Es evidente que a cualquier conjunto de celdas 
le corresponde cierto subconjunto de vectores y viceversa, 
a cualquier subconjunto de vectores le corresponde cierto conjunto 
de celdas. 

LEMA 1. Si el conjunto de vectores [P;,,) es linealmente depen- 
diente, el conjunto de celdas que le corresponde es ciclico. 

En efecto, supongamos que se sabe que el conjunto de vectores 
(P,,) es linealmente dependiente. Entonces existe una combinación 
lineal de estos vectores que lleva al vector cero q(0, 0, ..., 0) y por 
lo menos uno de los coeficientes es diferente de cero. Examinaremos 
sólo los vectores que entran en esta combinación con coeficientes 
diferentes de cero. Sea éste, por ejemplo, el vector P,,. Entonces 
en la combinación lineal obligatoriamente tiene que entrar por lo 
menos un vector con el mismo índice i,, por ejemplo el P,,. 
Eso tiene lugar porque en el vector q todos los componentes son 
iguales a cero, pero como en la combinación entra el vector P,,, 
para transformar el ¡-ésimo componente en cero hay que tener 
en la combinación aunque sea un vector más con el mismo 
componente diferente de cero. Puesto que ahora entra un vector 


ql 


que tiene el componente j,, entonces tiene que entrar por lo menos 
un vector más con el mismo componente, por ejemplo el vector 
P, ,  Razonando similarmente construimos la sucesión de los vectores 

O A O 


y? “ij? “lhajar A dapar 0 
Ya que se tiene sólo un número finito de vectores, después 
de cierto número tenemos que llegar al vector P,, y de él, al vector 


P, ;,- Esta sucesión de vectores corresponde al conjunto de celdas 
E) 1, 1172, l2j2, 123, qe.s Eje, 1, 
que es cíclico. 
LEMA 2. Cualquier conjunto compuesto de m + n celdas es ciclico. 
Para convencerse de la justeza de este lema es suficiente 
demostrar que cualquiera (m + n) vectores entre los construidos 
son linealmente dependientes. Demostrémoslo. 
Analicemos el vector (m + n)-dimensional. 


m componentes n componentes 


0 (—1, UL, ...> —1; +1, +1, ds + 1). 


Este es ortogonal con relación a todos los P;;. De aquí se 
deduce que todos los vectores P;; pertenecen al espacio (m + n— 
— 1)-dimensional y por eso cualquiera m + n de ellos son linealmente 
dependientes. 

De la condición (26”) del $ 11 y del último lema se deduce que 
n<N <m+mnm. En los ejemplos se puede advertir claramente que 
existen unos problemas para los que N =n y otros para los que 
N=m=>+mn- 1. Teniendo esto en cuenta, hacemos la conclusión 
definitiva de que N realmente satisface la condición n <N <m+ 
+n-=— 1 con lo que se demuestra el teorema. 

De los lemas que se acaban de demostrar se deduce que el 
conjunto de vectores (P;,,+ que corresponde a un conjunto acíclico 
de m+n-— 1 celdas es la base de los vectores P;, construidos 
en un espacio (m + n)-dimensional. Designaremos con la letra H 
un conjunto acíclico cualquiera de m+ n— 1 celdas. 

Teniendo en cuenta el mutuo y único significado de la relación 
entre todas las posibles bases del espacio (m + n)-dimensional 
y los conjuntos acíclicos de m + n— 1 celdas se pueden expresar 
las propiedades ya conocidas de las bases de los espacios 
multidimensionales aplicadas a los conjuntos acíclicos de m+n-— 1 
celdas de la siguiente manera: 

1” propiedad. Sea H, un cierto conjunto acíclico de m+mn-— 1 
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e (ij) H, *. Entonces el conjunto H,, obtenido al añadir (ij) al 
conjunto H,, contiene un ciclo O y sólo uno. 

2* propiedad. Sea (1, ¡YH (i, j) y (1, j) 8 O. Entonces el conjunto 
H3, que se obtiene de H, al excluir la celda (i', j'), de nuevo 
es un conjunto acíclico de m+n- 1. 

Los conjuntos H, y H3 sobre los que se trata en las 
propiedades 1 y 2 se diferencian solamente en una celda. 

Definición. Se llaman conjuntos acíclicos de m+mn-— 1 de una 
sustitución cada vez a dos conjuntos aciíclicos de m+mn-— 1 que 
se diferencian solamente en una celda. 

Examinemos una solución básica arbitraria X = (x;;). Según el 
teorema (véase la pág. 75) el número N de elementos x;; diferentes 
de cero satisface la condición n<N<m+mn-— 1. Los demás 
elementos m x n— N de la solución son iguales a cero. 

Construyamos un conjunto acíclico H de m+n-— 1 tal que 
todos los elementos x;; diferentes de cero de la solución X se 
encuentren en las celdas de este conjunto. 

Definición. Se llaman ceros seleccionados a los elementos x;; 
de la solución X iguales a cero y situados en las celdas del conjunto 
aciclico H de m+n- 1. 

Definición. Se llama selección al conjunto de elementos no cero 
xi; de la solución básica X junto con los ceros seleccionados que 
lo complementan hasta el conjunto acíclico H de m+n-— 1. 

Definición. Se llaman x-seleccionados a los elementos c;; de la 
matriz C que corresponden a los elementos de la elección X. 

En consecuencia, la selección es la solución básica en la que 
de todos los elementos cero x;¡=60 se han escogido aquellos que 
completan el conjunto de elementos no cero xj;X%0 hasta un 
conjunto acíclico H de m>+mn-— 1. Finalmente, llegamos a la 
deducción de que para encontrar aunque sea una solución óptima 
es suficiente investigar todas las posibles selecciones. 


$ 13. ELECCIÓN ÓPTIMA 


La metodología del cálculo de la solución óptima consiste 
en que, partiendo de cierta selección inicial, pasamos consecutiva- 
mente a otras selecciones con menor valor del producto escalar 
(CX) y después de un número finito de pasos llegamos a la 


*) El signo € significa que la celda ij no pertenece 
al conjunto H,. 
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solución óptima. Así surge la necesidad de formar la primera 
selección. 

Veamos el siguiente procedimiento para formar la primera 
selección: 

Determinamos los elementos de la primera línea de la matriz 
X = (x;;). Para ello buscamos y encontramos en la primera línea 
de la matriz C = (c;;) el elemento menor. Sea este el elemento Cc; . 
Entonces suponemos que x,,, = mín (a,;b;). Si a, > b; , encontramos 
en la misma línea el siguiente elemento menor que satisfaga 
a la condición C,; >C;, y suponemos x,, = mín (a; — x;;,; b;). 
Continuamos estos pasos hasta que no se satisfaga por completo 


n 


la primera ecuación a, = Y x,;. Si en cierto paso de este proceso 
La 


resulta que el resto de a, es exactamente igual al respectivo b, 


entonces tomando x,, igual a este resto x,, =b;,, suponemos 
completamente que el siguiente valor de la variable respectiva 
será x,;,,, =0. Después pasamos a la segunda línea, a la tercera, 
etc. En las columnas, cuyas ecuaciones están completamente 
satisfechas no se apunta el cero. 

De aquí se deduce que la primera solución, obtenida de esta 
manera, es una selección, puesto que el número de elementos 
marcados es igual a m+n-— 1 y entre ellos no hay ni una cadena 
cerrada. 

Veamos otro ejemplo para mayor claridad. Hay que formar 
la primera solución para un problema de transporte determinado 
con la matriz presentada en la tabla 3. 


Tabla 3 


En la primera línea, el menor de los números (8, 3, 5, 2) es el 
número 2, por eso tomamos x,Á = mín (10; 15) =10. Como la 
mercancía del primer punto de partida ya está distribuida por 
completo, pasamos a la distribución de la mercancía del segundo 
punto de partida. Para ello encontramos en la segunda línea 
de la tabla 3 el número menor (4, 1, 6, 7). El número 1 es el menor. 
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Por eso suponemos que x»» = mín (15; 10) = 10. Como en este caso 
queda resto en el segundo punto de partida igual a 15-10 =5 
unidades, buscamos en la segunda línea el siguiente elemento menor. 
Este será el número 4. Suponemos que x,, = mín (15 — 10; 5) =S5, 
Teniendo en cuenta que las necesidades del primer punto de destino 
están cubiertas por completo y que el resto en el segundo punto 
es igual a 0, encontramos el tercer elemento menor en la segunda 
línea. Tal elemento es el número 6. Por eso suponemos que x,3 =0. 
Ahora pasamos a la distribución de la mercancía del tercer punto 
de partida. Con este fin buscamos entre los números de la tercera 
línea (1, 9, 4, 3) el número menor. Como las necesidades del 
primer punto de destino están satisfechas por completo, no tomamos 
en cuenta en el cálculo el número 1 que se encuentra en el cruce 
de la tercera línea y la primera columna. Continuamos buscando 
el menor entre los números (9, 4, 3). Entre ellos el menor es el 
número 3. Suponemos que x34 = mín (25; 15 -— 10) = 5. Después 
encontramos en la tercera línea el siguiente número menor que 
es el 4 y suponemos que x33 = 20. 

Es fácil convencerse de que la distribución obtenida es la 
solución. Además la solución formada es selección. Efectivamente, 
el número de elementos elegidos es igualam+n-—1=3+4-1=6 
y ellos entre sí no forman ciclos. 

Se puede advertir fácilmente que, por el contrario, cualquier 
elemento no elegido x;; forma un ciclo y sólo uno con los elementos 
de la selección. 

A continuación emplearemos el procedimiento explicado para 
formar las primeras soluciones. 

Supongamos que se tiene la primera selección, la solución X, 
con el conjunto H,. Sea el respectivo producto escalar (CX;) 
igual a C,. Evaluemos cada elemento c;; que no entra en H, por 
medio de la expresión 


Aj= Ci LC 
o” o” 


en la que Y Ci, y Y Ci son las sumas de los elementos de las 
o" e» 
semicadenas impar y par, respectivamente, de la única cadena 
cerrada que forma cada elemento c;; con los x-seleccionados (el 
elemento evaluado c;; se toma como el primero en la cadena O). 
Señalemos el elemento c, , al que corresponde la menor evaluación 
Aín y su cadena cerrada con los elementos x-seleccionados. 
Formamos la segunda selección X, para lo que traspasamos el 
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elemento mínimo xp?" de la primera selección de la semicadena 
par a la semicadena. impar. En este caso, si en la semicadena par 
hay algunos elementos x;; iguales al menor xpn, entonces, para 
mayor certidumbre excluimos del conjunto H, la celda ij con el 
elemento xn que se encuentra la primera en la semicadena par 
al contornearla en la dirección de las agujas del reloj. En el lugar 
de la celda ij, excluida del as H,, introducimos la celda ¡,j, 
a la que le corresponde c,, y el nuevo elemento de la selección 
Xj,, = xppin (observemos que a cuenta del desplazamiento del 
elemento xmín por la cadena, en la selección X, pueden encontrarse 
varios elementos cero). Como la selección X, se diferencia de la 
selección X, solamente por los cambios en la cadena del elemento 
Cc, , , entonces es justa la desigualdad 


l1J1 
Y Ci, Loss 
O” 


el producto escalar (CX,) será menor que el (CX,) en el valor 
de 


DI gp y Ci] xp, 


o” o” 


Teniendo en cuenta que xpíh puede resultar igual a cero llegamos 
a la conclusión de que C, > C,. Se puede repetir todo el esquema 
de procedimientos realizado y partiendo de la selección X, formar 
la selección X3, etc. Como resultado se formará una sucesión 
de selecciónes X,, X», ..., Xy, .. tal que la respectiva sucesión 
de productos escalares C¡ >C7>..>C,>.. es una función no 
creciente del número de la selección. 

Se llama óptima la selección X, con el conjunto H, si cada 
uno de los elementos c,; que no entran en H, tiene un valor 
negativo A > 0. 

Puesto que para cualquier elemento c;,; que no entra en el 
conjunto H, es justa la desigualdad 


2 ij 2, Ci >0 
e” o” 


el paso, por medio de sustituciones únicas en cada vez, de la 
selección X, a cualquiéra otra selección no puede disminuir el valor 
del producto escalar con relación al producto escalar C; corres- 
pondiénte a la selección X, *. 


*) A continuación se demostrará que la selección 
óptima es al mismo tiempo la solución óptima. 
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TEOREMA. Para cualquier matriz C = (c;¡) de elementos reales 


en la que a; > 0, b;> 0, la sucesión de selecciones X,, X2, ..., Xy, ... 
después de un número finito de pasos se acaba con una selección 
óptima. 

DEMOSTRACIÓN. 

l-er caso. Si al pasar de una selección a otra, al formar 
la sucesión X,, X,, .. en ninguna de las selecciones aparece 


x¡¡¿ = 0, entonces xpín traspasado por la cadena también es diferente 
a cero. En este caso el valor del producto escalar en cada paso 
disminuye. Puesto que cuando m y n son finitos puede haber 
solamente un número finito de diferentes conjuntos H en cada uno 
de los cuales se define de una sola forma la selección, entonces 
no puede haber un número de pasos infinito. 

2-0 caso. Si al pasar de una selección a otra, en las nuevas pueden 
aparecer elementos x;, con valores cero, surge la incertidumbre 
sobre la calidad de finita de la sucesión X,, X,, .... Esta incertidumbre 
se basa en que en este caso xpiín=0 y el valor del producto 
escalar no disminuye. No obstante, esta incertidumbre se elimina 
por completo si simultáneamente al problema de transporte inicial 
se examina un problema con la misma matriz C = (c;;), pero con 


a; = 0; + €, E O 
b; = b;, Jj= 1, Ds 23ES (n E 1), 
b, = b, + me, 


donde e es un número positivo suficientemente pequeño. Llamemos 
a este segundo problema el e-problema. Formamos para estos dos 
problemas las primeras selecciones X, y Xí. Entonces como e 
es pequeñisimo, los conjuntos H, y Hí coincidirán y en la selección 
Xi no habrá elementos cero xj, =0. Efectivamente, la aparición 
de valores cero de x,¡ en una selección del problema inicial ocurre 
en los casos en los que existen combinaciones de líneas y columnas 
tales que se cumple la igualdad 


Por eso la aparición de valores cero x;, en el problema puede 
suceder en. dos casos: 


i=iq j=jp 
a) Na +qe= ) bj, 


i=1; j=j1 
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i=iq i=ip 
db) Y) a, +qe= Y b,+b, +me. 


i=1, =j1 


Escogiendo e diferente de las raíces de las ecuaciones expresadas 
en las condiciones a) y b) excluimos las posibilidades de que 
aparezcan elementos cero en la selección Xi. Dado que la matriz 
C =(c¡¡) es la misma para los dos problemas, coinciden los 
elementos C;¡, con la mayor evaluación negativa, las cadenas cerradas 
de estos elementos con los elementos de las selecciones X, y Xi, 
así como la situación en los conjuntos H, y Hi de los elementos 
xn y xgmín, destinados a ser transladados por la cadena. El paso 
de la selección X, a la selección X, en el problema inicial es 
de hecho el paso de la selección Xi a la selección X5 en el 
e-problema. Estos razonamientos, teniendo en cuenta que e es 
pequeñísimo, son correctos en cualquier paso. Como el g-problema 
es un problema en el que se carece de x;;=0 cuando tiene lugar 
la disminución del producto escalar en cada paso, pues en este 
e-problema después de un número finito de pasos llegaremos 
a la selección óptima X;. Debido a la coincidencia de las matrices 
C =(c;)) y los conjuntos H, y H;, la selección X, también es 
óptima. 
El teorema queda demostrado. 


$ 14. INVARIABILIDAD DE LA SUCESIÓN 
DE SELECCIONES, EQUIVALENTES A LAS TRANSFORMACIONES 
DE LA MATRIZ DE LOS COSTOS 


El cálculo de los elementos está ligado a un gran volumen 
de cómputo. Surge la tarea de definir transformaciones suficiente- 
mente sencillas de la matriz inicial C =(c;;), que sin cambiar 
la sucesión de selecciones, reduzcan considerablemente el proceso 
de búsqueda de la selección óptima. 

Para ello introduzcamos el concepto de transformación equi- 
valente de una matriz. 

Definición. Sea dada la matriz C = (c;¡) y los números arbitrarios 
Fiz Far. Fm 513 523 «3 Sn. Diremos que la matriz D = (d;) 
es equivalente a la matriz C = (c;;), si ella ha sido obtenida de 
la matriz C mediante la fórmula D = (c;;¡ + r; + sj). 

TEOREMA. La sucesión de selecciones es invariante con relación 
a las transformaciones equivalentes de la matriz. 
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DEMOSTRACIÓN. Transformemos la matriz C en una matriz D 
equivalente. Formamos para la matriz C la selección inicial Xy 


y partiendo de ella, por medio de una sola sustitución cada vez, 
formamos la sucesión X,, X,, X3, ... que converge hacia la 
selección Óptima X,. 

Partiendo de la misma selección X,, formamos para la matriz 
D la sucesión X,, X,, X3, .... convergente hacia la selección 
óptima X;. Fácilmente se puede ver que la diferencia de las 
evaluaciones A — AP? para cualesquiera dos elementos cf y cf 
que no entren en la selección X, es igual a cero: 


= [2 (cto — 20069] — 12 (es ¡E ries) 
puna 2 (có) — [2 (0) = 2 (ció) 


—Y (ri +5) +) (r,+5)=0 
O" e» 


Las últimas dos sumas son iguales entre sí y tienen diferente signo 
puesto que cualquier número r; (ó sj) contenido en la semicadena 
impar, se contiene también en la semicadena par. Como la 
evaluación de los elementos en las matrices C y D coincide, 
coinciden también los elementos que tienen el menor valor. Por 
eso, al pasar de la selección X, a la selección siguiente tanto 
para la matriz C como para la matriz D, las selecciones X, y X» 
coinciden. Estos razonamientos son justos en cualquier paso lo que 
quiere decir que la sucesión de selecciones X,, X», ..., Xp no 
cambiará al realizar transformaciones idénticas en la matriz. 
El teorema queda demostrado. 

Es evidente que esta afirmación no se altera al repetir las 
transformaciones idénticas en la matriz C. 

Teniendo en cuenta que los elementos de cualquier selección 
no forman cadenas cerradas entre sí, en cada paso se puede 
realizar una transformación idéntica que reduce a cero los elementos 
x-seleccionados. Para ello, por ejemplo, es suficiente añadir a los 
elementos de cada columna (línea) un número con el signo 
contrario al de suma algebraica del número que se añade a la 
línea (columna) y el valor del elemento x-seleccionado que se 
reduce a cero. 

La matriz con elementos cero x-seleccionados tiene la propiedad 
de que la evaluación de cualquiera de sus elementos que no sean 
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los x-seleccionados será igual al mismo elemento. Por eso no es 
necesario hacer grandes cálculos para determinar las evaluaciones 
de todos los elementos. Para formar la selección siguiente es 
suficiente marcar el elemento negativo mayor (según su valor 
absoluto). Si ya todos los elementos de la matriz transformada 
resultan no negativos al mismo tiempo que los elementos 
x-seleccionados son cero, entonces la última selección es ya la 
Óptima. De la invariabilidad de la sucesión de las matrices se 
deduce que la selección óptima es al mismo tiempo la solución 
óptima. Efectivamente, la solución óptima con la matriz trans- 
formada de manera equivalente, en la que los elementos 
x-seleccionados son iguales a cero, tiene la propiedad de que su 
producto escalar es menor (o por lo menos igual) que el producto 
escalar de cualquiera de las soluciones. Por esa razón la selección 
Óptima que hemos encontrado es la solución óptima. 


$ 15. ALGORITMO DEL CÁLCULO DE LA SOLUCIÓN ÓPTIMA 


I. Los datos del problema se anotan en tablas. 

Il. Se define la primera selección. 

III. Se reducen los elementos x-seleccionados a cero. Si después 
de eso los demás elementos de la tabla transformada resultan 
no negativos, entonces la primera selección ya es la óptima. 

IV. Si después de la reducción a cero de los elementos 
x-seleccionados en la tabla hay números negativos, encontramos 
el elemento negativo de mayor valor absoluto. 

V. Formamos una cadena cerrada del elemento negativo mayor 
con los elementos x-seleccionados reducidos a cero y componemos 
la segunda selección transpasando por la cadena el número x; 
el que sea el menor en la semicadena par. 

VI. Como el elemento negativo mayor ha resultado un 
x-seleccionado, lo reducimos a cero, pero de tal manera que queden 
iguales a cero los demás elementos x-seleccionados que corres- 
ponden a la nueva selección. 

VIT. Repetimos estos pasos hasta que lleguemos a una selección 
tal en la que en la tabla transformada todos los elementos 
x-seleccionados sean iguales a cero y los demás sean no negativos. 
Esta última selección será la solución óptima. 

VII. Para calcular los gastos de transporte que corresponde 
a la solución óptima, se deben reunir en una misma tabla la matriz 
inicial de costos y la última solución. La suma de los productos 


J> 
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de los dos números que se encuentran en cada una de las celdas 
de la tabla determina el valor de los gastos de transporte. 

Para mayor claridad, se presenta la interpretación geométrica 
del proceso de cálculo de la solución óptima (fig. 20). 


C(A,) - costo inicial 


A, C/(A,) - costo final 


0 A Ar Aj Ay As As 4, 


Fig. 20 


Marquemos en el eje horizontal los números de las soluciones 
y en el eje vertical, los respectivos valores de los gastos. 

Al comenzar, formamos la primera solución y encontramos uno 
de los puntos, el A,. La reducción a cero de los elementos 
x-seleccionados nos lleva al punto A' del eje horizontal. Si en 
este punto la tabla transformada no tiene elementos negativos, 
la solución que corresponde al punto A, es la óptima. Si la tabla 
transformada no tiene esta propiedad, encontramos el punto A»); 
de nuevo volvemos al eje en el punto 4% y comprobamos si hemos 
alcanzado ya la solución óptima. Si no la hemos logrado todavía, 
encontramos el punto Ax, etc. 

Después de un número finito de pasos se llegará forzosamente 
al último punto, el que corresponde a la solución óptima. 
Obtendremos la línea quebrada de los gastos que aparece en la 
fig. 20. Como se puede ver en el dibujo, después de cada paso 
el valor del gasto disminuye (o en todo caso no aumenta). 

Continuemos la resolución del problema de transporte planteado 
en la pág. 79. En la tabla 4 se presenta la primera selección para 
este problema. Reduzcamos a cero los elementos x-selecciónados. 
Para ello restemos 4 de los elementos de la primera columna 
de la matriz de costos; de los elementos de la segunda columna 
restemos 1. Como resultado de esto se reducirán a cero los 
elementos x-seleccionados que se encuentran en los cruces de la 
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segunda línea con la primera y la segunda columna. Restemos 6 
de los elementos de la tercera columna, entonces el elemento 
x-seleccionado que se encuentra en el cruce de la segunda linea 
con la tercera columna también se reducirá a cero y el elemento 
de la matriz de costos que se encuentra en la intersección de 
la tercera línea con la tercera columna tomará un valor igual 
a (4-6= -—2). Añadamos a los elementos de la tercera línea 
de la matriz de costos el número 2 y de los elementos de la 
cuarta columna restemos el número 5. Entonces, los elementos que 
se encuentran en los cruces de la tercera línea con la tercera 
y la cuarta columnas se reducirán a cero. 


Tabla 4 


20 


pa 
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Ea 


Tabla 5 


Por último, queda reducir a cero el elemento x-seleccionado 
que se encuentra en el cruce de la primera línea con la cuarta 
columna. Para ello hay que añadir a los elementos de la primera 
línea el número 3. Como resultado de esta transformación idéntica 
de la matriz de costos, la tabla 4 se convertirá en la tabla S. 

Como se ve en la tabla 5, hay un elemento negativo (— 1). 
Por eso no podemos estar seguros de que la primera selección 
es la óptima. Como en este ejemplo hay un solo elemento negativo, 
el mismo es el de mayor valor absoluto. 

Formemos una cadena cerrada de este elemento con los 
x-seleccionados y compongamos la segunda selección. Eso se 
muestra en la tabla 6. Con este fin transpasemos de la semicadena 
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par a la impar el número de unidades de mercancías que en la 
semicadena par sea el menor. 

Como el número —1 ya es x-seleccionado, hay que reducirlo 
a cero, para lo cual sometemos de nuevo la matriz de costos 
a transformaciones idénticas. Añadimos a los elementos de la 
primera columna de la tabla 6 un número igual a 1. Obtendremos 
como resultado la tabla 7 en la que no hay ni un elemento 
negativo, y todos los elementos x-seleccionados son iguales a cero. 


Tabla 6 


La selección presentada en la tabla 7 es, basándonos en lo 
indicado, la solución óptima. Reunimos en una sola tabla 8 
la matriz de costos inicial y la solución óptima. El gasto de 
transporte que corresponde a este plan es C=2-10+5-6+ 
+ 10-145-344.15+1-5= 140 unidades de costo. 

Según este plan hay que organizar el transporte de la 
siguiente manera (fig. 21): Toda la mercancía del primer punto 
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de partida se transporta al cuarto consumidor; el segundo punto 
de partida abastece con 10 unidades de mercancía al segundo 
consumidor y con 5, al tercero; desde el tercer punto de partida 
se tansportan 5 unidades de mercancías al primer consumidor, 
15 al tercero y las 5 unidades de mercancías restantes se envían 
al cuarto consumidor. Con este plam, el gasto de transporte 
(a igualdad de las demás condiciones) será el mínimo. 


Obras N4 
O) 15 


Fábrica ws 


10 MÍ 


LA 
Fabrica n2 
15 

O) 
Obras N?7 
Í 


Fig. 21 


Examinemos un ejemplo más complicado. 

Ejemplo. Hallar el plan de transporte óptimo de una mercancía 
homogénea desde cuatro puntos de partida hasta seis puntos 
de destino para las condiciones definidas en la tabla 9. En esta 
tabla a, y b, se miden en miles de toneladas y c;;, en miles 
de rublos. 


Tabla 9 
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Escribamos otra vez la tabla inicial, introduzcamos en la parte 
Inferior de las celdas los números de la selección inicial y presentemos 
todo el proceso de obtención de la solución óptima como una 
sucesión de las tablas 9a, b, c, d, e, f. 


Tabla Ya 


Reducimos 
Primera a cero los 
selección elementos 
C,=75 x—seleccionados 
de la tabla 
Tabla 9b 
E 
Se marca con Traspasamos 
un cuadrado (0) 5 unidades por 
el elemento que la cadena. 
se introduce da El gasto dis— 
la solución O : 
dd minuye en 
al 20 unidades 
O 
AN 
+4 +4 +4 
Tabla 9c 
El elemento que se excluye de la 
solución aparece sombreado 
Segunda Transpasamos 
selección O unidades por 
la cadena. 


C,=55 El gasto queda 


el mismo 
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Tabla 9d 


Primera cadena 


ALMAS 


Transpasamos por 
la primera cadena 2, 


Tercera por la segunda 1. 
selección DS (0/0) El gasto disminuye 
C,=55 ÍA en 2x2+2x1=6 


unidades 


a TOS 


+2 =£ 


Segunda e 


Tabla Ye 
Cuarta Transpasamos una 
selección unidad es cadena. 
C=4 El gasto disminuye 
es en 1x1=1 unidad 
Quinta Solución 
selección Optima 
C¿=48 


Unamos la tabla inicial de costos 
y la solución Optima (tabla 10). 


En la última tabla se indica cómo hay que organizar el 
transporte, para garantizar el mínimo de gastos. La variante 
óptima en comparación con la primera variante garantiza una 
disminución de los gastos de transporte de 27 mil rublos, ya que 
estos gastos en la primera variante son de 75 mil rublos y en 
la óptima sólo de 48 mil rublos. 
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Tabla 10 


Es de interés observar que la tabla que contiene la variante 
óptima (la 5? selección) se indica no sólo la única solución 
Óptima, sino todas las soluciones con los mismos gastos. 

Efectivamente, en esta tabla hay elementos cero que no están 
marcados con un círculo. Marquémoslos con una estrella. Cada 
uno de estos elementos forma una cadena cerrada con los elementos 
de la solución. 

Empleando con las reglas anteriores el procedimiento de una 
sola sustitución cada vez se pueden obtener otras soluciones que 
tienen los mismos gastos. 

Esto permite tomar como óptima una solución, por ejemplo, 
como la presentada en la tabla 11. Los gastos correspondientes 
a este plan de transporte también son iguales a C=1-2+4-3 + 
+4-14+1-54+1:-14+ 1-24 1-24+5.1+4+3.5=48. 

Fácilmente se ve que en este caso cada cadena cerrada compuesta 
sólo de elementos cero, tiene la propiedad de que son iguales las 
sumas de los gastos en las semicadenas par e impar de la 
respectiva cedena cerrada de la tabla inicial lo que permite 
traspasar xpín por la cadena sin alterar el gasto total. En una 
serie de casos resulta muy útil definir todas las soluciones óptimas. 
Como veremos más adelante, aplicando al conjunto de soluciones 
óptimas el procedimiento de elección de la variante óptima de 


Tabla 11 


93 


transporte bajo el criterio del tiempo, se puede encontrar un plan 
de transporte que se realice en un tiempo mínimo. 

Por último, refirámonos a la composición de la cadena cerrada 
formada por el elemento negativo mayor con los elementos 
x-seleccionados. Cuando m y n no son muy grandes, las cadenas 
se ven directamente como en los ejemplos examinados anterior- 
mente. En cambio, si m y n son relativamente grandes, la búsqueda 
de la cadena cerrada se realiza de la forma más sencilla de la 
manera siguiente: se tachan las líneas y las columnas que contienen 
un cero marcado (puesto que para componer la cadena hay que 
tener no menos de dos ceros en la línea o en la columna); 
después de tachar la primera, se tacha la segunda, etc., hasta que 
lleguemos a una situación en la que enseguida se vea la cadena 
cerrada. 

Por cuanto la resolución detallada de un ejemplo que tenga 
una tabla de grandes dimensiones nos ocuparía mucho lugar, aquí 
nos limitaremos a examinar un caso en el que se debe hallar la 
cadena cerrada. Esto se presenta en la tabla 12. Al principio 
tachamos todas las columnas que contienen un cero (claro está, 
excluyendo la columna en la que se encuentra el elemento negativo 
mayor). Después examinemos todas las líneas y tachamos aquellas 
en las que ha quedado un cero. De nuevo revisamos las columnas 
y hacemos lo mismo. 

Ahora para componer la cadena cerrada solamente hay que 
unir los elementos que quedan comenzando por el número negativo 


Tabla 12 
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y revisando los elementos cero, por ejemplo, en la dirección 
de las agujas del reloj. 

Hemos estudiado la metodología de resolución de los problemas 
de transporte bajo el criterio del costo a condición de que 


m n 
Y a; = » b,. En este caso, la cantidad de mercancía en todos 
¡=1 j=1 

los puntos de partida es igual a la cantidad de mercancía necesaria 


en todos los puntos de destino. 
S1 


entonces introducimos un punto de destino ficticio con un consumo 
m n 
br + LA y ad; — y b 
i=1 j=1 


y con un costo de transporte en este punto Cc; ,-¡=0 para 
O A 

Ejemplo. En cuatro depósitos hay respectivamente 70, 40, 50 
y 20 toneladas de combustible. Se debe planificar el transporte 
de combustible a tres consumidores de tal manera que el gasto 
del transporte sea mínimo. En la tabla 13 se dan las condiciones 
del problema. En este caso es preferible formar la primera 
solución por las columnas. Esta solución está marcada con 
círculos en la tabla 14. 

Es fácil comprobar que la solución presentada en la tabla 14 
es la óptima. De acuerdo con esta solución al primer consumidor 
hay que llevarle 10 toneladas del primer depósito y 20 toneladas 
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del cuarto; al segundo consumidor, 30 toneladas del primer 
depósito y 10 toneladas del segundo; al tercer consumidor, 
30 toneladas del segundo depósito; además en el primer depósito 
quedan aún 30 toneladas de combustible; en estas condiciones 
no es conveniente emplear las reservas del tercer depósito. 


Composición de la selección inicial 


Reducción a cero de los elementos x-se- 
leccionados de la matriz de los costos 
(transformación de la matriz) 


Búsqueda en la matriz transformada del 
elemento negativo de mayor valor absoluto 


Control del final de la resolución: se com- 
prueba la existencia de elementos negati- 
vos en la matriz transformada | 


Cálculo del costo de 
transporte según el plan 
óptimo 


Salida a la impresora 


Composición de la cadena cerrada formada 
por el elemento negativo de mayor valor 
absoluto y los elementos x-seleccionados 


Composición de una nueva selección 


Fig. 22 


del plan de transporte 
óptimo y el valor del 
costo correspondiente 


El algoritmo expuesto para la resolución del problema de 
transporte bajo el criterio del costo puede presentarse como un 
proceso rutinario de operaciones aritméticas y lógicas. Este proceso 
se realiza fácilmente en una máquina computadora electrónica. 
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En la fig. 22 se ilustra un esquema de bloques para la resolución 
de este problema. A base de este esquema de bloques fue preparado 
un programa para la resolución del problema en la computadora 
“Strelá”. 

Para una matriz de un orden de mxn=<500 el tiempo de 
resolución del problema no es mayor de 8—10 minutos; para una 
matriz de un orden de m x n < 1500 pueden hacer falta para resolver 
un problema hasta 25—40 minutos; para una matriz de un orden 
m x n de varios miles hacen falta hasta 8—10 horas para resolver 
un problema en la computadora “Strelá”. 

Por ejemplo, en la computadora “Strelá” se resolvió el siguiente 
problema de transporte. 

Desde nueve puntos de partida hay que transportar una 
mercancía a catorce puntos de destino (tabla 15). La selección 
inicial obtenida en la computadora está representada en la tabla 16. 


Tabla 15 
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Tabla 16 
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Tabla 17 


/ 2 3 4 5 8 10 "12 E 


e OODOoDonEa 


EX 
EN 
aaa 
BEN 
NA 


e 
E ¡E 
EN 


7 


Y 13) 
(e) e 
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Era 
E jor 
EA ES 
ae 55 
A Ea 
ES —— E 
E MINA 
E AAA 
EA A 04 10 0000 1500 0 E A 


La selección final que corresponde a la solución óptima se ilustra 
en la tabla 17. 

Es fácil calcular que los gastos que corresponden a la selección 
inicial son de un valor igual a C =971 y los que corresponden 
a la óptima, C = 703, 

El tiempo de resolución de este ejemplo (en un solo cálculo) 
fue de menos de un minuto. 

Otro problema de transporte con una matriz de costos de 
mx n=30 x 38 (solucionada en dos cálculos) fue resuelto en un 
lapso de 37 minutos. 


E A A 
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CAPÍTULO IV 


RESOLUCIÓN DEL PROBLEMA DE TRANSPORTE 
POR EL CRITERIO DEL TIEMPO 


La economía de tiempo, al elaborar un plan de transporte, es 
absolutamente necesaria en una serie de casos prácticos importantes. 
Por ejemplo, al transportar productos de fácil deterioro es necesario 
su traslado a los puntos de destino en el plazo mínimo posible. 
En la época de la cosecha un factor muy importante es el rápido 
traslado del grano a los centros de acopiamiento. 

Los problemas de este tipo son los de transporte por el criterio 
del tiempo. A continuación se estudia uno de los algoritmos de 
la resolución del problema de transporte por el criterio del tiempo. 


$ 16. PLANTEAMIENTO Y RESOLUCIÓN DEL PROBLEMA 


Sean dados m puntos de partida de una mercancía homogénea 
y n puntos de destino. Sean a;,, 4, ...., Qi, ... Am la cantidad 
de unidades de mercancía que se encuentran en al primero, segun- 
do, ..., i-ésimo, ..., m-Ésimo puntos de partida y b,, b,, ..., b;, ..., by 
la cantidad de unidades de mercancía que tienen que ser entregadas 
en cada uno de los n puntos de destino. Sea t;; el tiempo 
(en días o en horas) necesario para el transporte de la mercancía 
desde el i-ésimo punto de partida hasta el j-ésimo punto de destino. 
Sea x;¡¡ la cantidad de unidades de mercancía que planificamos 
transladar del ¡-ésimo punto de partida al j-ésimo punto de 
destino. Hay que hallar el plan de transporte óptimo, es decir, 
unos números no negativos x,; con los que el tiempo de transporte 
de todas las mercancías necesarias a todos los puntos de destino 
sea el mínimo. 

El planteamiento matemático de este problema se reduce 
a lo siguiente: 

Se da un sistema de ecuaciones algebraicas lineales: 


(27) 


99 


en las que los números a, y b, satisfacen la condición 
m n 
24= ) Dj 
¡=1 j=1 


Sea, además, dada la matriz del tiempo T = (t;;). A cada solución 
no negativa X =(x;) del sistema de ecuaciones (27), es decir, 
a cada plan de transporte le corresponde una determinada matriz 
Tx = (tj), cuyos elementos son tf = t;; si x¡¿>0 y tf =0 si x;¡¿ =0. 
Llamemos té el elemento mayor de la matriz Tx. 

Se debe hallar una solución no negativa X¿p,, para la cual 
o sea el menor entre todos los £** que corresponden a diferentes 
soluciones no negativas X. Es cómodo presentar las condiciones 
de este problema en forma de la tabla 18. 


Tabla 18 


Como se puede mostrar fácilmente, este problema no se resuelve 
con la ayuda del algoritmo que fue descrito para el problema 
de transporte por el criterio de gastos. 

Por ejemplo, la solución para la que la forma lineal C = ) t;;x;; 
alcanza el valor mínimo, como regla general, no es la óptima por 
el tiempo. Para el ejemplo dado en la tabla 19 e ilustrado en 
la fig. 21 se presenta la solución óptima por el mínimo de la 
forma lineal C = Y t;,x,¡ en la tabla 20. En la tabla se ve que toda 
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la mercancía será trasladada hasta los puntos de destino solamente 
después de seis días. Por otra parte, la variante de transportes 
presentada en la tabla 21 corresponde a un cierto valor mayor 
de la forma lineal, pero el traslado de las mercancías a los puntos 
de destino se cumple en cuatro días. 

Al transportar productos de fácil deterioro los gastos en el 
recorrido adicional de vagones (camiones) se cubren con la 
conservación de la calidad de miles de toneladas de productos 
alimenticios destinados para la población. 

La demostración del método de resolución que se describe más 
adelante se deduce de la teoría general de la programación lineal 
que fue expuesta en el capítulo Il. 

Veamos la descripción general del método en la resolución del 
problema que se plantea con la tabla 22. En esta tabla se ve que 
hay que transportar 125 unidades de mercancías que se encuentran 
en seis puntos de partida a siete puntos de destino. Supongamos 
que los números de esta tabla significan el tiempo necesario, 
en horas, para el transporte de la carga desde el respectivo punto 
de partida hasta los de destino. Por ejemplo, el número 31 que 
se encuentra en el cruce de la cuarta línea con la cuarta columna 


Tabla 20 Tabla 21 
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muestra que el tiempo de transporte de la carga desde el cuarto 
punto de partida hasta el cuarto punto de destino es 31 horas. 
Hay que hallar el plan de transporte óptimo por el tiempo. 

Representemos la tabla anotada (llamémosla la tabla pequeña) 
en forma de una tabla grande (tabla 23). 

La línea superior contiene los índices de las variables x;; desde 11 
hasta 67, la inferior, los valores de t;; contenidos en las respectivas 
celdas ij. Por ejemplo, el par de números 43 que están en la línea 
de arriba significa el lugar de la variable x¿3, el valor del tiempo 
taz que corresponde a esta variable es 39; por eso 43 y 39 se 
encuentran en una columna. 

Toda la tabla está dividida en zonas. El número de zonas es 
igual al número de puntos de partida o, lo que es lo mismo, 
al número de lineas en la tabla inicial. Además, la tabla está 
dividida horizontalmente. En la mitad superior de las zonas hay 
seis líneas: el número de líneas en la tabla inicial; en la parte 
inferior, hay siete líneas: el número de sus columnas. En la 
semizona superior al nivel de cada número figuran siete signos 
menos en cada respectiva zona vertical. Por ejemplo, al nivel del 
30, los menos están anotados en la cuarta zona. Si examinamos 
la tabla inicial veremos que esta anotación significa la expresión 
del resto de la carga en el cuarto punto de partida: 


SÓ — Maz — Xa2 — Xaz — Xag — Xas — Xa6 — Xa7- 


Similarmente a esto, en la parte inferior la colocación de los signos 
menos en cada línea expresa la falta de mercancía en los puntos 
de destino. Por ejemplo, la colocación de los menos en la línea 
donde está el 27 indica que la falta de mercancía en el punto 
de destino que corresponde a la cuarta columna es igual a 


lO O SA CE PP O CE 


Con otras palabras, la colocación de los menos en la tabla 
corresponde a las ecuaciones que expresan las condiciones del 
problema. | 

Construyamos la primera solución utilizando el procedimiento 
empleado en el capítulo III. Distribuyamos las 15 unidades del 
primer punto de partida entre los puntos de destino hasta los que 
el tiempo del recorrido es el menor; las 7 unidades del segundo 
punto de partida, entre los puntos de destino que todavía no 
están abastecidos hasta los que el tiempo del recorrido es el menor, 
etc., etc. Asi se forma la primera solución presentada en la 
tabla 24 que, en general, puede no ser la óptima. Como se ve, 
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Tabla 24 
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esta composición elemental nos ha llevado a una solución cuya 
realización exige 28 horas para el transporte de la mercancía. 
Eso nos permite, a continuación, excluir de la revisión las celdas 
que tienen t;¡> 28, lo que se ha hecho en la tabla anexa. 

Para la composición de la siguiente variante de solución, cuya 
realización exija menos de 28 horas, dirijámonos a la tabla grande. 
Para los datos de esta tabla, la composición de la selección se 
hace de la siguiente manera: en la línea inferior de los valores t;; 
se halla el menor de los números de la primera zona. Ese es el 
número 7 que corresponde a la columna 14, o sea, a la variable x;4. 
Examinando esta columna nos convencemos de que podemos hacer 
X14 Igual a 15, ya que la reserva de mercancía es de 15 unidades 
y la necesidad, de 27. Eso permite anotar la variable x,¿ a la 
izquierda de la igualdad en la primera línea y suprimir el signo 
menos que se encuentra en el cruce de la primera línea y la 
columna 14. El segundo menos de esta columna que se encuentra 
al nivel del número 27, también se puede suprimir, y de la línea 
que corresponde al número 27 se resta la línea que corresponde: 
al número 16. 

En las líneas primera y décima de la primera zona han habido 
cambios. Estas líneas tienen ahora la siguiente forma: 


x14 =15/1—|—|—1|-]-]|=|. 


= 12 |+|+|+]|+]|+]+1. 


Pasamos a la segunda linea y a la segunda zona. Hallamos 
el menor de los números t;; de la segunda zona. Ese es el número 
6 que corresponde a la variable x2¿. Revisando la columna 
verticalmente, elegimos el menor de los dos números que se 
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encuentran a la izquierda al nivel de los signos menos de esta 
columna. Ese es el número 5. Anotemos la variable x,6 a la 
izquierda de la igualdad enfrente del 5 suprimiendo el menos que 
se encuentra en el cruce de esta línea y la columna 26. Suprimimos 
también el menos del cruce de la segunda línea y la columna 26, 
y de la segunda línea que corresponde al número 7 restamos 
la línea correspondiente al número 5. En la segunda línea queda 
el número 2. Eso significa que del segundo punto de partida el 
producto no está todavía distribuido entre los puntos de destino. 
Por eso, en la línea de los valores de t;; de la segunda zona 
hallamos el número siguiente por su valor. Ese es el número 7 
que se encuentra en la columna de la variable x,,. Revisando 
la columna encontramos el menor de los números de la izquierda 
que se encuentran al nivel de los menos de la columna (20, 2). 
Suprimimos los menos en la columna 21, anotamos la variable x,, 
a la izquierda de la igualdad enfrente del 2 y restamos de la 
séptima línea la segunda. Con esto acaba la distribución de las 
unidades de la segunda línea. Continuamos este proceso hasta 
la distribución de las unidades de la sexta línea. Como resultado 
obtendremos la primera solución que coincide con la solución 
anotada en la tabla 24. La tabla grande se ha transformado 
ahora en la tabla 25. 

De la justeza de la solución nos convencemos por la falta 
de los “+” y los “—” en una de las líneas. En esta línea, 
en el caso de que no haya errores en la distribución, tiene que 
haber sólo ceros. En nuestro ejemplo, la última línea ha resultado 


m U 


ser de este tipo. En virtud de la condición ) a,= ) bj, una de 
¡=1 j=1 

las ecuaciones del sistema (27) debe ser consecuencia de las demás. 

Las columnas que corresponden a las variables básicas en la 
tabla están señaladas con un rayado vertical y las columnas en las 
que su f;; > 28, con un rayado inclinado. Estas últimas se excluyen 
de la revisión posterior. En la tabla 25 están indicados los modos 
de composición de una mejor variante en comparación con la 
primera. Realmente, nosotros deseamos hallar una solución tal en 
la que la variable que se encuentre en la misma columna que 
t;¡¡= 28, es decir la x;7, sea igual a cero. Ahora bien, el contenido 
de la sexta línea de esta tabla: muestra que la disminución del 
valor de la variable x¿, puede ser obtenida a cuenta del aumento 
de xs, o de xs. Por cuanto a la variable x¿, le corresponde 
te, =15 y a la variable xós, tos = 23 haremos la disminución 
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de xs, a costa del aumento de x¿,. Revisando la columna 61 
vemos que el menor número que se encuentra en la parte izquierda 
de la tabla enfrente del “—” de esta columna es el valor xg, = 5. 

Ahora realicemos la transformación de la última tabla de la 
siguiente manera: excluyamos las columnas que tienen el rayado 
inclinado y la columna que corresponde a la variable x¿, (o bien, 
lo que es lo mismo, al número t;;, igual a 28). En lugar de xs, 
escribamos x¿,. Siguiendo, de las líneas que tienen signo menos 
en la columna 61 restamos la línea que corresponde a la variable 
Xe61, y a las líneas que tienen signo más en la columna 61, les 
sumamos la línea que corresponde a esa misma variable. Después 


de esto desaparecerán todos los “*—” y los “+” de la columna 61. 
Tabla 26 

12/14/1517 21[22/26/27|31|52|34/35 41/42[45/67|51| 57/61 
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Como resultado obtendremos la tabla 26. De esta tabla se 
deduce que ya hemos compuesto la segunda solución para la 
realización de la cual se necesita bastante menos tiempo: en lugar 
de 28 horas, sólo 21. (Observemos que si hubiésemos comenzado 
la reducción de la variable x¿, a cero no a costa del aumento 
de xs,, sino del aumento de xs, hubiese sido compuesta una 
variante de solución con un tiempo máximo = ts = 23. Después, 
tachando todas las columnas con t;;> 23, hubiésemos tenido que 
hacer otro paso para componer una tercera solución en la que 
el máximo de tiempo fuese t;¡= 21). 

Enunciemos las variables, el cambio de cuyo valor ha llevado 
a la segunda solución: X67, X35, X47, X31, Xas, Xó61- Estas variables 
han formado una cadena cerrada (tabla 27). Numeremos las celdas 
con las variables que forman la cadena tomando como primera 
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la celda con la variable x¿, (la que se tiene que excluir). Esta 
cadena tiene la propiedad de que en sus celdas impares infalible- 
mente se encuentran elementos x-seleccionados y en las pares, los 
valores  t;¡¡<ti?%, donde ti” es el valor de t;¡¿ que debe de ser 
excluido (en este caso t2f'! = 28). Llamemos aligeradora a esta 
cadena. Después de hallar la cadena aligeradora para formar una 
nueva solución traspasamos por ella una cantidad de unidades 
de carga que sea igual a la mínima entre ellas en la semicadena 
impar. En la tabla pequeña esta composición nos lleva a una 
segunda solución que coincide con la segunda solución compuesta 
en la tabla grande. Excluyamos de la investigación posterior las 
columnas (celdas) con valores t;¡> 21. Así obtendremos la tabla 
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pequeña 28 y la tabla grande 26 (contando que las columnas 41, 
55 y 65 están tachadas). Surge una pregunta: ¿habrá alguna solución 
que se realice en menos de 21 horas? Para dar contestación 
a esta pregunta examinamos la tabla 26. En la línea de los 
valores t;; hallamos el mayor de ellos, el 21. A este valor t34 = 21 
le corresponde la variable x3¿, en la línea de la cual hay sólo 
signos “+”. La falta de los “—” indica la imposibilidad de 
disminuir el valor de x34 a cuenta de alguna otra variable. 
No se ven posibilidades de excluir x34 y, por consecuencia, 
de liberarse de tz, = 21. Así resulta que la última variante obtenida 
es la óptima por el criterio de tiempo. Como se dijo anteriomente 
esta variante se realiza en 21 horas. En este caso no existe 
la posibilidad de componer una cadena aligeradora en la tabla 
pequeña 28. 

Para aclarar mejor el método de cálculo de la variante óptima 
por el criterio del tiempo, hemos empleado tablas grandes 
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y pequeñas. En realidad, al resolver tales problemas hace falta 
sólo una tabla grande o pequeña hecha a lápiz. Los pasos de una 
variante a otra se efectúan en una tabla empleando lápiz y goma 
de borrar. 


Tabla 28 


C=1716K 


Las reglas para el cálculo de la variante óptima por el criterio 
del tiempo cuando se emplea una tabla pequeña se reducen 
a lo siguiente: 

1. Anotamos las condiciones del problema en una tabla pequeña. 

2. Encontramos la primera solución (por ejemplo, con el proce- 
dimiento indicado anteriormente, 

3. Determinamos el t¡¡máx que corresponde a esta solución. 

4. Tachamos en la tabla todas las celdas que tengan t;; > t¿¿máx, 

5. Investigamos la posibilidad de componer cadenas aligeradoras 
formadas de t¡¡máx con los elementos que quedan en la tabla. 

6. Si aparecen cadenas aligeradoras componemos una solución 
nueva. 

7. Si no es posible aligerar por completo la celda que 
corresponde a t¡¡máx (reducir la respectiva x,, a cero) entonces 
la solución que contiene t¡¡náx es la óptima. 

8. Si se aligera por completo la celda que contiene t¡¿máx, 
entonces hallamos t¿máx en la nueva solución obtenida. 

Así hemos cumplido el primer paso dirigido a la obtención 
de la solución óptima. Si en este caso no se ha llegado a la 
solución Óptima entonces hay que hacer un segundo paso, 
cumpliendo de nuevo las reglas anotadas desde el punto 4. 

Para obtener la solución óptima habrá que realizar un número 
de pasos finito. 

Al emplear estas reglas, la única dificultad consiste en la 
composición de las cadenas aligeradoras. Con respecto a esto, 
el empleo de la tabla grande, cuyas reglas de transformación 
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fueron descritas anteriormente con suficiente detalle, da la 
posibilidad de llegar sin grandes impedimentos a la solución 
óptima. 

Se podrían exponer también otros procedimientos de resolución 
de este problema. 

Cuando m y n son pequeños, estos problemas pueden ser 
resueltos con papel y lápiz. Si m y n son grandes entonces es 
necesario, lo mismo que en los problemas que se resuelven por 
el criterio del costo, el empleo de las máquinas computadoras 
electrónicas. En este caso el tiempo de máquina de resolución 
de estos problemas, para los mismos valores de m y n, coincide 
aproximadamente con el tiempo de resolución de problemas por 
el criterio del costo. 


$ 17. RESOLUCIÓN DE LOS PROBLEMAS 
DE TRANSPORTE TOMANDO EN CUENTA 
EL TIEMPO Y EL COSTO 


En la práctica se presentan condiciones en las que es más 
conveniente la realización de un plan de transporte que sea un 
promedio entre el óptimo por tiempo y el óptimo por costo. 
Supongamos, para hacer más sencilla la descripción, que el costo 
es proporcional al tiempo de transporte, o sea C¡;= Kt;;. 

En la tabla 29 se presenta la variante óptima, por el criterio 
de costos, del problema visto anteriormente que se cumple en 
28 horas. Los gastos de transporte de esta variante son C= 
=K (7-154+6-5+10-2+11.-20+20-13+4+21-124+5.9+ 12-21+ 
+ 14.12 + 16-11 + 28-5) = K- 1668 unidades de costo. Los gastos 
en el transporte que corresponden a la variante óptima por el 
criterio del tiempo (tabla 28) son de C=K(7:-154+7-2+6-5+ 
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+ 11-13 + 20-13 + 21-12 + 21-74 5-2 + 12-28 + 14-12 + 
+ 15.5+ 16-11) = 1716K unidades de costo. Haciendo un gasto 
complementario, igual a 1716K — 1668K = 48K acabaremos las 
operaciones de transporte de mercancías 7 horas antes en 
comparación con el plan de transporte determinado por la variante 
óptima, calculada por el criterio del costo. 

Del ejemplo examinado se deduce que, admitiendo un gasto 
complementario en el transporte relativamente pequeño, se puede 
reducir considerablemente el tiempo de las operaciones de transporte 
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de los productos. En los casos en los que se dan los plazos del 
cumplimiento de las operaciones de transporte se puede, empleando 
los métodos descritos, hallar la variante óptima por el criterio 
de costos que se realice en el margen del plazo dado. Si la pérdida 
de tiempo es inadmisible, entonces calculamos la variante óptima 
sólo por el criterio del tiempo. Si la variante obtenida con tgPtno 
es la mejor por los costos, entonces empleando el algoritmo del 
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cálculo de la solución Óptima por el mínimo de la forma lineal, 
podremos obtener la óptima por los costos entre todas las óptimas 
por el tiempo (o sea, entre todas las que se realizan en un mismo 
plazo, igual a tg). Como se ve en el ejemplo examinado, la 
solución óptima por el tiempo presentada en la tabla 30 no es 
óptima por los costos. Efectivamente, reduzcamos a cero los 
elementos x-seleccionados (al reducir a cero los elementos 
x-seleccionados, no se toma en cuenta el valor de los elementos 


Tabla 32 


que se encuentran en las celdas tachadas, tabla 31) y como resultado 
obtendremos una tabla en la que hay un número negativo (—14). 
Transformando otra vez la tabla 31 obtendremos la tabla 32. 
En la tabla 32 ya transformada los elementos x-seleccionados son 
iguales a O y los demás son no negativos. Por eso, la variante 
de solución compuesta es la óptima por el criterio del costo. 
El costo de la realización de esta variante, que se lleva a cabo 
también en 21 horas, es de C =1688K lo que es mucho menos 
que 1716K. 


* *K * 


En los últimos años, los métodos de programación lineal 
encuentran cada vez mayor aplicación en la economía, la técnica, 
el arte militar, y en otras esferas. La teoría y los métodos de 
programación lineal se perfeccionan constantemente y permiten 
solucionar otros nuevos problemas. El desarrollo acelerado de la 
técnica de computación ha hecho posible la solución práctica de 
cualquier problema de programación lineal. 

La continuación del desarrollo de los métodos de programación 
lineal y su empleo para la resolución de problemas de la economía 
nacional va a favorecer la mejora de la organización y de la 
planificación de la industria nacional. 
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